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An  asymptotic  expansion  method  for  a nonlinear  singularly  perturbed  two-point 
boundary  value  problem  arising  in  a nonlinear  optimal  control  problem  with  free  final 
time  is  presented.  Assuming  that  i)  the  reduced  solution  satisfies  Jacobi’s  sufficient 
condition  for  a local  minimum  for  the  reduced  problem  and  ii)  the  fast  subproblem 
satisfies  an  eigenvalue  criterion,  the  asymptotic  validity  for  the  asymptotic  solution  is 
proved  via  a Banach  space  fixed  point  theorem.  The  proposed  asymptotic  expansion 
method  is  extended  to  a nonlinear  two-point  boundary  value  problem  with  multiple 
small  parameters. 

As  a design  example,  the  proposed  method  is  applied  to  the  derivation  of  a 
guidance  law  for  an  air-to-air  bank-to-turn  missile.  In  order  for  the  asymptotic 
expansion  to  remain  valid  during  the  entire  interval,  a control  weighting  adjustment 


scheme  is  used.  The  resulting  guidance  law  is  easily  implementable  for  an  on-line 
computation  and  gives  a good  performance  under  both  the  deterministic  and  the 
stochastic  environments. 
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CHAPTER  I 
INTRODUCTION 


The  development  of  optimal  control  theory  in  the  early  1960s  brought  promises 
of  many  new  areas  for  applications.  In  order  for  a solution  of  an  optimal  control 
problem  to  be  optimal,  the  solution  must  satisfy  the  necessary  condition  expressed 
in  a two-point  boundary  value  problem  (TPBVP).  But  a TPBVP  resulting  from 
application  of  Pontryagin’s  minimum  principle  to  an  optimal  control  problem  is  in 
general  very  difficult  to  solve  analytically  except  some  special  cases,  for  example, 
a linear  system  with  a quadratic  cost.  Accordingly,  many  numerical  methods  have 
been  developed  to  solve  a nonlinear  TPBVP  effectively  [1]  [2] . 

However,  in  a certain  class  of  systems  with  small  parameters  multiplying  deriva- 
tives, the  resulting  TPBVP  is  stiff  so  that  a numerical  method  may  fail  to  solve  the 
TPBVP.  This  type  of  system  is  singularly  perturbed  and  frequently  appears  in  con- 
trol problem  formulations.  For  example,  when  one  considers  i)  a dynamic  system 
with  small  parameters  such  as  small  masses  and  small  inertias  etc.,  ii)  a dynamic 
system  with  a fast  actuator  and  iii)  a high  gain  feedback  system,  the  dynamic  system 
can  be  formulated  as  a singularly  perturbed  system. 

One  intuitive  and  widely  used  approach  to  design  such  a control  system  is  to 
neglect  small  parameters  or  fast  dynamics.  The  basic  reason  for  this  approach  is 
that  one  can  design  a control  system  from  a lower  order  model.  But  in  this  case, 
the  closed  loop  system  designed  from  the  lower  order  model  may  be  unstable  or 
the  performance  of  the  control  system  may  be  unacceptable  [3].  Another  approach 
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is  to  seek  a trade-off  between  the  lower  order  system  and  the  high  order  system 
involving  small  parameters  or  fast  dynamics  such  that  the  design  by  a trade-off  gives 
a better  performance  than  that  designed  from  the  lower  order  model  and  takes  less 
computational  burden  comparing  with  the  design  from  the  high  order  model.  Many 
authors  have  used  the  singular  perturbation  method  to  achieve  this  goal  during  the 
past  two  decades  [4]. 

The  singular  perturbation  method  has  been  a widely  used  technique  to  obtain 
an  approximate  solution  in  the  area  of  an  applied  mathematics  such  as  fluid  me- 
chanics and  celestial  mechanics.  The  singular  perturbation  method  has  been  mainly 
applied  to  an  optimal  control  problem  since  the  singular  perturbation  method  was 
first  introduced  in  a control  problem  formulation  by  Kokotovic  [5].  But,  little  has 
been  done  for  a nonlinear  optimal  control  problem  with  free  final  time.  In  a special 
class  of  nonlinear  system  whose  fast  state  is  scalar,  some  efforts[6]  - [9]  was  made 
to  derive  a zeroth  order  or  a first  order  solution  from  an  optimal  control  with  free 
final  time  formulation,  but  neither  an  explicit  method  to  construct  an  asymptotic 
solution  nor  its  asymptotic  validity  was  treated. 

In  this  dissertation,  we  have  two  objectives.  One  is  to  develop  an  explicit 
method  to  construct  an  asymptotic  solution  in  a nonlinear  optimal  control  problem 
with  free  final  time  and  fixed  end  points  which  plays  an  important  role  in  a design 
of  an  aircraft  guidance  law.  The  other  is  to  derive  a guidance  law  for  an  air-to-air 
bank-to-turn  (BTT)  missile. 

The  major  control  issues  in  designing  a BTT  missile  guidance  law  are  i)  the  BTT 
missile  motion  equation  is  nonlinear,  ii)  the  guidance  law  must  be  implementable  by 
on-board  computation,  iii)  the  allowable  control  is  limited  to  a control  constraint  and 
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iv)  the  only  target  information  available  are  the  relative  distances  in  each  axis  and 
the  closing  velocity  corrupted  by  noise.  The  guidance  problem  may  be  formulated 
as  an  optimal  control  problem.  Since  the  missile  motion  equation  is  nonlinear,  the 
resulting  nonlinear  TPBVP  is  difficult  to  solve  by  on-line  computation. 

The  utility  of  the  singular  perturbation  method  when  applied  to  a nonlinear 
system  is  that  an  approximate  solution  can  be  obtained  by  solving  several  lower 
dimensional  subproblems  instead  of  solving  the  full  problem.  Hence,  in  order  to 
successfully  apply  the  singular  perturbation,  the  problem  must  be  formulated  such 
that  the  lower  dimensional  subproblem  is  solvable  without  heavy  computational 
burden.  In  this  respect,  recently,  the  singular  perturbation  method  has  been  succes- 
sively utilized  in  obtaining  a missile  guidance  law  [6] [7].  But,  it  was  observed  that 
the  nearly  optimal  guidance  law  derived  based  on  the  singular  perturbation  method 
failed  to  give  a good  approximation  during  the  terminal  phase  in  which  the  time 
scale  separation  is  not  valid  [7]. 

In  the  application  of  optimal  guidance  laws  to  missile  control,  feedback  guid- 
ance laws  must  be  used  to  compensate  for  target  motion,  measurement  errors  and 
model  errors.  However,  this  means  that  the  optimization  problem  should  be  resolved 
many  times  as  the  time-to-go  (final  time  minus  initial  time)  changes.  This  causes  a 
difficulty  with  the  singular  perturbation  scheme  in  that  the  choice  of  a performance 
index  which  allows  time  scale  separation  depends  on  time-to-go.  Therefore,  a per- 
formance index  which  will  allow  good  time  scale  separation  and  corresponding  good 
approximations  for  large  time-to-go  will  fail  to  work  for  small  time-to-go. 

In  this  dissertation,  a simple  method  to  specifically  overcome  the  difficulties 
encountered  in  the  application  of  the  singular  perturbation  technique  to  generating 
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optimal  feedback  guidance  laws  over  varying  time  intervals  is  developed.  The  actual 
goal  of  the  guidance  law  is  to  minimize  the  miss  distance,  but  the  performance  index 
which  is  introduced  also  includes  a weighting  on  the  control  effort  to  avoid  impulsive 
controls.  The  control  weighting  influences  the  time  scale  separation  of  the  solution 
which  is  an  important  assumption  in  the  singular  perturbation  technique.  A method 
for  adjusting  the  control  weighting  of  the  performance  index  will  be  presented  such 
that  a uniform  degree  of  approximation  is  obtained  as  the  time-to-go  approaches 
zero.  This  modification  leads  to  a significant  improvement  in  the  final  miss  distance. 

The  outline  of  this  dissertation  will  be  as  follows. 

In  Chapter  II,  the  problem  is  set  up  and  notation  necessary  for  the  development 
for  a singularly  perturbed  free  final  time  optimal  control  problem  is  introduced  and 
previous  results  on  singular  perturbation  theory  are  reviewed. 

In  Chapter  III,  an  explicit  method  to  construct  an  asymptotic  solution  is  de- 
scribed. An  asymptotic  solution  will  be  constructed  in  three  different  time  scales. 

In  Chapter  IV,  it  will  be  shown  that  the  asymptotic  solution  constructed  in 
Chapter  III  is  asymptotically  valid  via  a Banach  space  fixed  point  theorem. 

In  Chapter  V , a nonlinear  system  with  multiple  small  parameters  will  be  con- 
sidered. For  a nonlinear  system  with  multiple  parameters,  we  will  construct  an 
asymptotic  solution  by  extending  the  method  developed  in  Chapter  III  to  a multi- 
ple parameter  case. 

In  Chapter  VI,  the  main  drawback  of  the  currently  widely  used  guidance  law 
will  be  reviewed  and  by  using  the  method  described  in  Chapter  V,  a BTT  guidance 
law  up  to  the  first  order  will  be  constructed. 
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In  Chapter  VII,  the  simulation  results  will  be  illustrated  for  the  case  when  i) 
all  the  measurements  are  available  under  non-ideal  autopilot/air  frame  dynamics 
(deterministic  case)  and  ii)  only  the  bearing  only  measurements  are  available  under 
non-ideal  autopilot/air  frame  dynamics  (stochastic  case). 

The  final  chapter  gives  conclusions  and  some  suggestions  for  further  works. 


CHAPTER  II 

BACKGROUND  AND  PROBLEM  STATEMENT 


2.1  Brief  Summary  of  Singular  Perturbation  Theory 


In  this  section,  a mathematical  tool  developed  in  order  to  obtain  an  approx- 
imate solution  for  a differential  equation  with  a multiple  time  scale  property  is 
reviewed.  Singularly  perturbed  systems  and,  more  generally,  multiple  time  scale 
systems,  often  occur  due  to  the  presence  of  small  parasitic  parameters  such  as  small 
time  constants,  masses  and  capacitors  etc.  The  major  advantages  of  the  singular 
perturbation  approach  in  an  analysis  of  a system  are  that  one  can  obtain  a good  ap- 
proximate solution  from  a lower  order  system  and  avoid  an  ill  conditioning  problem 
resulting  from  stiffness  of  the  system  [5]. 

First,  some  preliminary  notation  will  be  defined.  We  are  interested  in  the  limit 
of  functions  such  as  /(e)  as  e goes  to  zero.  Instead  of  saying  that  /(e)  goes  to  zero 
at  the  same  rate  that  e goes  to  zero,  we  say  /(e)  is  order  e or  write  it  as  /(e)  = 0(e). 
In  general,  we  put 


/(e)  = 0(sr(e))  as  e ->  0 if  lim4r  = * 

t-o  g(e) 

for  some  finite  constant  k.  Let  (5n(e)  be  a general  sequence  of  functions.  Sn(e)  is 
called  an  asymptotic  sequence  if  it  satisfies 


lim 

e— >0 


<Wi(e) 

Sn(e) 


= 0. 
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Given  Y^=oan^n{^)i  where  the  an  does  not  depend  on  c and  <5n(e)  is  an  asymp- 
totic sequence,  we  say  that  this  expansion  is  an  asymptotic  expansion  and  write 

OO 

/(c)  ~ £an5„(c). 

n=0 

We  remark  that,  given  an  asymptotic  sequence  <5n(e),  the  representation  of  /(e) 
in  terms  of  this  sequence  is  unique  [10].  In  our  development,  en  will  be  used  as  a 

<5»(e). 

We  will  define  |x|,  a norm  of  a vector  or  a matrix  x,  as  the  sum  of  the  absolute 
values  of  all  the  elements  of  x.  Given  f{x(t),y(t),z(t)),  fx(t ) will  be  used  to  de- 
note the  partial  derivative  of  f(x(t),y(t),z(t))  with  respect  to  x(t)  unless  otherwise 
specified.  Thus, 

fx(t)  = | £{x(t),y(t),z(t)). 

The  superscript  T denotes  the  transpose  of  a vector  or  a matrix. 

For  convenience,  the  following  notational  convention  will  be  used.  If  x denotes 
any  vector  in  1Z2n  for  some  integer  n,  then 


x - 


x\ 

x2 


where  xx  and  x2  denote  the  first  n and  the  last  n components  of  x,  respectively. 
Also,  for  2 n x 2 n matrix  A,  denotes  n x n submatrix  for  1 < i,j  < 2.  Thus, 


A = 


An  A12 
A21  A22 

First  of  all,  we  will  consider  a simpler  problem,  namely  the  singularly  perturbed, 


uncontrolled,  autonomous,  initial  value  problem  which  was  treated  by  numerous 
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authors  [11]  - [13]. 

x = 
ez  = 

with  the  initial  condition 

x(0,e)  = x0 


f(x,z) 

(2.1) 

g{x,z) 

(2.2) 

O 

II 

'V 

CD 

where  x(t,e)  and  z(t,e)  are  n and  m dimensional  vectors,  respectively,  and  e is  a 
small  positive  constant  which  represents  the  ratio  of  the  time  derivative  of  x and  z. 
The  initial  conditions  Xq  and  Zq  are  n and  m dimensional  constant  vectors. 

We  are  interested  in  (2.1)  - (2.2)  under  the  assumption  that  e is  small  relative 
to  the  other  parameters  of  the  system.  Since  i is  large  compared  with  x in  this  case, 
we  will  call  x the  slow  state  and  2 the  fast  state.  It  is  reasonable  to  set  e = 0 when 
we  want  to  obtain  an  approximate  solution.  Setting  e to  zero  gives 


f{xr,zr) 

(2.3) 

d(Xr,Zr) 

(2.4) 

with  the  initial  condition  xr(0)  = x0. 

We  will  call  the  system  (2.1)  - (2.2)  the  full  problem  and  the  system  (2.3)  - 
(2.4)  the  reduced  problem  associated  with  the  full  problem.  It  can  be  observed  that 
the  order  of  the  reduced  problem  is  n.  Thus,  by  perturbing  e to  zero,  the  order  of 
the  system  is  changed  from  n + m to  n.  For  this  reason,  this  type  of  perturbation 
is  called  singular.  But  in  the  reduced  problem  (2.3)  - (2.4),  in  general,  zr  does  not 
satisfy  the  initial  condition  since  zr  is  simply  a solution  of  the  algebraic  equation 
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(2.4).  In  order  to  circumvent  this  difficulty,  the  time  scale  is  stretched  by  introducing 
a stretched  time  scale  r = ^ and  the  boundary  layer  terms  x L(r)  and  zl(t).  In  the 
stretched  time  scale,  the  full  system  will  be 

dxL 


dr 

dzL 

dr 


= e/(A  z1) 
= »(z£',zi) 


(2.5) 

(2.6) 


with  the  initial  conditions 


xL(0,e)  = xo  zL(0,e)  = zo. 

The  system  (2.5)  - (2.6)  is  called  the  boundary  layer  problem  and  the  time 
interval  r E [0,  oo)  is  called  the  boundary  layer.  In  order  to  obtain  an  approximate 
solution  of  the  problem,  it  is  natural  to  set  e = 0 in  (2.5)  - (2.6).  Then,  the  resulting 
system  becomes 


dzL 

dr 


= x0 

= g(x0,zL) 


(2.7) 

(2.8) 


with  the  initial  condition  zL( 0)  = z0. 

It  can  be  observed  that  the  solution  of  the  reduced  problem  at  t = 0 is  an 
equilibrium  point  of  the  system  (2.7)  - (2.8).  One  crucial  condition  for  the  validity 
of  this  approximation  scheme  is  the  stability  of  the  system  (2.7)  - (2.8).  We  will 
state  a theorem  taken  from  Kokotovic’s  work  [14]. 

Theorem  2.1.  Suppose  that 


1-  f{xiz)  and  g(x,z)  are  continuously  differentiable  in  some  open  region  Q of 


their  domain. 
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2.  The  reduced  problem  has  a unique  solution  on  some  interval  0 < t < T. 

3.  There  exists  an  isolated  root  zr  = <f>(xT)  of  (2.4)  in  Q. 

4.  The  solution  zl(t)  of  (2.8)  exists  on  r £ [0,  oo)  and  is  unique.  Moreover,  the 
system  (2.8)  is  asymptotically  stable  with  respect  to  the  root  <f>(xr)  of  (2.4). 

5.  The  real  parts  of  the  eigenvalues  of  the  Jacobian  matrix 

§9,  \ 

dz^Xr,Zr^ 

are  negative  on  [0,  T]  for  zT  = <f>(xr) 

Then,  for  sufficiently  small  e,  the  full  system  (2.1)  and  (2.2)  has  a unique  solution  on 
t 6 [OjT]  satisfying  the  initial  conditions  i(0,c)  = xo  and  z(0,e)  = zq.  Furthermore, 

limx(f,  e)  = xr(t) 

uniformly  on  0 < t < T 

limz(«,e)  = zr(t) 

uniformly  on  any  closed  subinterval  of  0 < t < T. 

If  a more  accurate  approximate  solution  than  is  obtained  in  the  reduced  prob- 
lem is  required,  an  asymptotic  expansion  method  can  be  used.  This  method  was 
pioneered  by  Vasil’eva  [15].  The  basic  idea  is  that  a solution  is  sought  in  the  form 
of  a sum  of  asymptotic  expansions  which  approximates  the  solution  both  inside  and 
outside  the  boundary  layer. 

For  various  kinds  of  TPBVPs,  an  asymptotic  expansion  method  and  its  va- 
lidity was  investigated  by  many  authors  [16]  - [22].  Kokotovic  and  Sannuti  first 
applied  the  singular  perturbation  method  to  a control  system  design  in  order  to 
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obtain  a reasonably  good  approximate  solution  from  a lower  order  model  [5].  After 
that,  Wilde  and  Kokotovic  also  applied  a singular  perturbation  method  to  a lin- 
ear regulator  problem  [16].  The  necessary  conditions  for  regulator  problem  for  a 
singularly  perturbed  linear  system  leads  to  a singularly  perturbed  Riccati  equation 
which  also  has  the  multiple  time  scale  property.  Therefore,  an  approximate  solution 
can  be  obtained  by  solving  a lower  order  Riccati  equation  and  an  algebraic  equa- 
tion. O’Malley  also  applied  the  singular  perturbation  method  to  a linear  regulator 
problem  with  a quadratic  performance  index  [17].  He  then  applied  the  singular 
perturbation  method  to  the  resulting  linear  TPBVP. 

Hadlock  considered  a class  of  nonlinear  TPBVP  arising  in  an  optimal  control 
problem  with  a fixed  final  time  and  free  end  points  [23].  lie  showed  that  the  slow 
states  and  costates  of  the  full  problem  converge  to  the  slow  states  and  costates  of  the 
reduced  problem  as  e goes  to  zero,  uniformly  on  the  entire  closed  interval,  whereas 
the  fast  states  and  costates  of  the  full  problem  converge  to  the  fast  states  and  costates 
of  the  reduced  problem  on  the  open  interval  (0,T).  The  lack  of  convergence  of  the 
fast  state  and  costate  at  the  end  points  results  from  the  fact  that  the  fast  states  and 
the  costates  of  the  reduced  problem  which  are  the  solution  of  an  algebraic  equation 
do  not  satisfy  the  boundary  conditions.  But  Hadlock  did  not  give  an  explicit  method 
to  construct  an  asymptotic  power  series  solution. 

For  a so-called  quasi-linear  system  in  which  the  slow  states  are  nonlinear  but 
the  fast  states  and  controls  are  linear,  O’Malley  developed  a procedure  to  construct 
an  asymptotic  power  series  solution  for  an  optimal  regulator  problem  [18].  Based  on 
Hadlock’s  work,  Sannuti  developed  an  explicit  method  to  construct  an  asymptotic 
power  series  solution  and  showed  asymptotic  validity  via  a successive  approximation 
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scheme  in  a class  of  nonlinear  TPBVPs  [19].  But  Sannuti  also  restricted  the  TPBVP 
to  a problem  resulting  from  a fixed  final  time  and  free  end  points  quasi-linear  optimal 
control  problem.  Given  a quasi-linear  optimal  control  problem  with  a fixed  final 
time  and  a fixed  end  point  problem,  Chow  developed  a method  to  construct  a 
power  series  solution  [20].  He  constructed  the  boundary  layer  terms  by  solving  an 
algebraic  Riccati  equation. 

For  a general  nonlinear  TPBVP  resulting  from  a nonlinear  fixed  final  time  free 
end  point  optimal  control  problem,  Freedman  developed  a procedure  to  construct  an 
asymptotic  series  solution  uniformly  valid  throughout  the  interval  [0,T]  and  proved 
its  asymptotic  validity  [21]  [22].  Recently,  Visser  and  Shinar  developed  a procedure 
to  construct  an  approximate  solution  in  feedback  form  in  a special  class  of  nonlinear 
optimal  control  problems  with  free  final  time  [9].  They  assumed  that  i)  a solution 
from  the  reduced  problem  can  be  obtained  in  feedback  form,  ii)  the  fast  state  is 
scalar,  iii)  there  exists  no  boundary  layer  at  the  terminal  point  and  iv)  the  Hamilto- 
nian is  not  an  explicit  function  of  time.  The  main  motivation  for  their  work  is  that 
one  can  obtain  a boundary  layer  solution  by  using  the  fact  that  the  Hamiltonian  is 
identically  zero.  In  this  case,  one  does  not  have  to  solve  a TPBVP  in  the  boundary 
layer.  But  they  did  not  prove  the  validity  of  their  method. 

Breakwell  and  Shinar  developed  an  alternative  method  to  construct  a feedback 
solution  in  a class  of  nonlinear  optimal  control  problem  with  free  final  time  [24]. 
They  used  the  Hamilton- Jacobi-Bellman  equation  instead  of  using  the  resulting 
TPBVP.  Their  method  will  be  briefly  reviewed. 

Consider  the  system 


x = f(x,z,u ) x(0)=aro 


(2.9) 
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ei  = g(x,z,u)  2(0)  = z0  (2.10) 

where  x is  an  n dimensional  slow  state,  z is  a scalar  fast  state  and  u is  an  r dimen- 
sional control  variable. 

The  objective  is  to  bring  the  system  from  the  given  initial  condition  to  a ter- 
minal manifold  determined  by  x(T ) while  minimizing  the  performance  index 

J — / L(x,z,u)dt 
Jo 

where  T is  an  unspecified  final  time. 

Define  the  optimal  return  function  V(x(t),  z(t),  e) 

, , rT 

V(x(t),z(t),e)  = min  J | L(x,z,u)dt' . 

Then,  along  an  optimal  trajectory,  the  Hamilton-Jacobi-Bellman  partial  differ- 
ential equation  must  be  satisfied 

0 = min  (L(x,z,u)  + Vj  f{x,z,u)  + ^-g(x,z,u)).  (2.11) 

The  optimal  control  u will  be  a function  of  x,z,  Vx  and  \VZ. 

u = U(x,z,Vx,-Vz). 

e 

By  setting  e = 0 in  (2.9)  and  (2.10),  (2.9)  and  (2.10)  become 

x — f(x,z,u) 

0 = 9(x,z,u). 

Let  2 = <f>{x,u)  be  the  solution  of  0 = g(x,z,u).  Then  the  reduced  problem  is 
given  by 

x = f(x,<t>(x,u),u)  (2.12) 


0 = g(x,</>(x,u),u). 


(2.13) 


14 


Assume  that 

1.  The  optimal  return  function 

[T 

Vr0(x(<))  = min  J L(x,<f>(x,u),u)dt 
can  be  obtained  for  the  reduced  problem. 

2.  The  optimal  return  function  in  the  full  problem  has  an  asymptotic  expansion. 
Thus, 

OO 

V(x(t),z(t),e)  ~ VJel 

3=0 

Substituting  the  asymptotic  expansion  V(x(t),  z(t),  e)  into  (2.11)  and  matching 
the  coefficients  of  equal  powers  of  e will  give,  for  k > 1, 

1/  ^ n 1 

Vk,z  = = hk(x,z) 

where  hk(x,z ) is  a function  of  x(<),z(t),  VjtX  and  VhZ  for  0 < j < k — 1. 

Since  hk(x,z)  is  a known  function  and  z(t)  is  a scalar,  Vk  can  be  obtained  by 
integrating  hk(x,z ) with  respect  to  z(t).  Hence  one  can  obtain  correction  terms  in 
the  optimal  return  function  recursively  as  long  as  integration  can  be  done  analyti- 
cally. Even  if  this  method  requires  the  optimal  cost  function  of  the  reduced  problem 
and  analytical  integration,  there  are  many  applied  problems  such  that  this  method 
can  be  used. 


2.2  Problem  Statement 

In  this  section,  we  will  consider  a nonlinear  singularly  perturbed  TPBVP  arising 
in  an  optimal  control  problem  with  free  final  time.  The  problem  formulation  is 
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similar  to  the  problem  studied  by  Breakwell  and  Sliinar  [24];  however,  we  will  present 
a different  analytical  approach.  Consider  the  following  nonlinear  system. 


X = f(x,z,u ) 

x(0)  = x0 

(2.14) 

tz  = F(x,z,u) 

u 

0 

II 

£ 

(2.15) 

where  x and  z are  n and  m dimensional  state  vectors,  respectively,  u is  an  r dimen- 
sional control  variable  and  e is  a small  positive  constant. 

The  objective  is  to  obtain  a control  u(t,e)  which  minimizes  the  following  scalar 
performance  index  with  the  terminal  constraint. 

J = + L(x,z,u)dt  (2.16) 

x{T)  = xT  (2.17) 

where  T is  free,  xT  is  prescribed,  L(x,z,u)  is  nonnegative  and  limr^  <f>(T)  = oo. 

The  functions  /,  F and  L can  be  functions  of  e.  However,  it  will  be  assumed 
that  /,  F and  L are  not  functions  of  e for  simplicity.  In  order  to  obtain  the  necessary 
conditions,  we  will  introduce  the  Hamiltonian 

H(x,  z,  p,  q,  u)  = L(x,  z,  u)  + pTf(x,  z , u)  + -JT F(x , 2,  u)  (2.18) 

where  p and  q are  the  costate  variables  corresponding  to  the  state  variables  x and 
2. 

For  convenience  of  singular  perturbation  analysis,  it  is  customary  to  use  q = -q 
as  the  fast  costate  variables  instead  of  using  q [17]  — [22].  Application  of  Pontryagin’s 
minimum  principle  gives  the  following  nonlinear  TPBVP  which  must  be  satisfied 
along  an  optimal  trajectory. 


P = ~ Hx(x,z,p,q,u ) 


(2.19) 


e<7  = ~ Hz(x,z,p,q,u ) q(T)  = 0 


(2.20) 


0 = Hu(x,z,p,q,u) 

(2.21) 

0 = [<M0  + H(x,z,p,q,u)]t_T. 

(2.22) 

Now,  suppose  that  u = U{x,z,p,q)  satisfies  the  equation  (2.21).  Substitution 

°f  u = Z'P'  *?)  'nto  (2-14)  - (2.15)  and  (2.19)  - (2.22)  gives  the  following  2n  + 2m 

dimensional  TPBVP. 

x = f(x,z,p,q)  z(0)  = x0  x(T)  = xT 

(2.23) 

p = g{x,z,p,q) 

(2.24) 

ez  = F(x,z,p,q ) z(0)  = z0 

(2.25) 

eq  = G(x,z,p,q)  q(T)  = 0 

(2.26) 

0 = [<M0  + H(x,z,p,  q)]t=T 

(2.27) 

where 

f(x,z,p,q)  = f(x,z,p,q,U(x,z,p,q)) 

(2.28) 

9(x,z,p,q)  = -Hx(x,z,p,q,U(x,z,p,q)) 

(2.29) 

F(x,z,p,q)  = F{x,z,p,q,U(x,z,p,q )) 

(2.30) 

G{x,z,p,q)  = -Hz(x,z,p,q,U{x,z,p,q)) 

(2.31) 

H(x,z,p,q ) = H(x,z,p,q,U(x,z,p,q)). 

(2.32) 

By  setting  t — 0,  the  reduced  problem  can  be  obtained  as  follows. 


xr  = f(xr,zr,pr,qr) 


(2.33) 
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Pr 

S{xri  zriPrtQr) 

(2.34) 

0 

F(xr,Zr,Pr,<]r) 

(2.35) 

0 

— G(xr,zr,pr,qr) 

(2.36) 

0 

+ H(xr,  zr,pr,  (lr)]t—xr 

(2.37) 

*r(0) 

Xr(Tr)  — X'j'. 

(2.38) 

One  can  observe  that  the  reduced  problem  consists  of  a 2 n dimensional  differ- 
ential equation  and  a 2m  dimensional  algebraic  equation.  In  computing  the  solution 
of  the  full  problem  by  a numerical  method,  the  computational  burden  will  be  pro- 
hibitively large  as  the  order  of  the  full  problem  increases.  Since  the  full  problem 
is  stiff,  a small  integration  step  size  is  required  in  order  for  the  integration  to  be 
numerically  stable,  further  increasing  the  computational  burden.  But  in  the  reduced 
problem,  these  difficulties  may  be  eliminated  and  it  is  possible  to  obtain  an  analytic 
solution  for  a special  class  of  nonlinear  TPBVPs.  The  most  important  questions  in 
the  study  of  singular  perturbation  are 

1.  Under  which  conditions  will  the  solution  of  the  full  problem  approach  the 
solution  of  the  reduced  problem  as  e goes  to  zero? 

2.  Under  which  conditions  can  we  obtain  a uniformly  valid  approximate  solution 
on  the  interval  [0,T]  from  the  reduced  problem? 

3.  How  can  we  make  correction  terms  in  order  to  obtain  a more  accurate 
approximate  solution? 

The  above  questions  for  the  class  of  nonlinear  TPBVPs  described  in  this  section 
will  be  investigated  in  this  dissertation. 


CHAPTER  III 

ASYMPTOTIC  EXPANSIONS 


In  this  chapter,  we  will  construct  an  asymptotic  solution  of  the  full  problem 
(2.23)  - (2.27)  based  on  the  reduced  solution.  The  construction  will  be  performed 
in  three  different  time  scales.  The  construction  method  is  similar  to  the  earlier 
methods  [13]  [19]  [21]  for  a fixed  final  time  problem  except  that  the  final  time  also 
can  be  expanded  into  a power  series  in  e. 


3.1  Basic  Approach 

Since  the  solutions  of  the  reduced  problem  do  not  satisfy  the  initial  and  the 
terminal  conditions,  we  can  expect  that  the  solution  of  the  full  problem  will  have 
initial  and  terminal  boundary  layer  terms.  Also,  we  can  expect  that  the  optimal 
final  time  will  be  a function  of  a small  parameter  e.  Moreover,  we  will  assume  that 
all  the  state,  costate  variables  and  the  optimal  time  can  be  expanded  into  a power 
series  in  e. 

It  will  be  assumed  that  all  the  variables  can  be  expressed  in  the  following  form. 


x(t,e)  = x°(t,e)  + cxl(t,  e)  + exR(a,e) 

(3.1) 

p(t,e)  = p°(t,e)  + epL(r,e)  + epR(a,e) 

(3.2) 

z(t,e)  = z°(t,e)  + zL(T,e)  + zR{<r,e) 

(3.3) 

q(t,e)  = q°(t,e)  + qL(T,e)  + qR(a,e) 

(3.4) 
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where 
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t T -t 

t = -,  <r  = 

e e 


We  seek  a solution  by  means  of  the  following  asymptotic  power  series. 


( x°(i,£),P°(t,t),z0(t,e),q0(t,e ))  ~ 

OO 

E(x?  (o»  z°  (o»  <■. i ? wy 

3=0 

(3.5) 

( xL(t , e),pL(r,  e),  zl(t,  e),  qL(r,  e))  ~ 

OO 

izocj'(T)’Pj  (r)»  2jl(t)’  9j'(T))ej 

3=0 

(3.6) 

{xR(a,  e),pR(a,  e),  z*(cr,  e),  gH(cr,  e))  ~ 

OO 

E(xf  (a)’Pf  (<t)>  *?(*)»  (cr))eJ 

3=0 

(3.7) 

OO 

r~  EV 

3=0 

(3.8) 

3.2  Outer  Expansion 


We  will  seek  a solution  (x°(t,  e),p°(t,  e),  z°(t,  e),  q°(t,  e))  which  represents  the 
behaviour  of  the  solution  of  the  full  problem  in  the  region  away  from  the  end  points. 

Assumptions. 

AJ^  f(x,  z,p,  q),g(x,  z,p,  q),F(x,  z,p,  q),G(x,  z,p,  q),  H(x,z,p,q ) and  <£(T)  de- 
fined in  (2.28)  - (2.32)  are  infinitely  differentiable  with  respect  to  all  their  arguments. 

A.2.  The  reduced  problem  (2.33)  - (2.38)  has  a unique  and  infinitely  differen- 
tiable solution  on  the  interval  0 < t < Tr. 

A.3.  The  reduced  solution  satisfies  the  sufficient  conditions  [25]  for  a local  min- 
imum. Thus, 


1.  The  convexity  condition 
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2.  The  normality  condition 

3.  The  no  conjugate  point  condition. 


AA  The  2m  x 2m  matrix  B(t)  is  invertible  for  all  t £ [0,Tr],  where  B(t)  is  defined 

by 


F, (t)  Fq{t) 

G, {t)  Gq(t) 


and,  for  example, 

dF 

Fz(t)  = -fc(xr(t),zT(t),pr(t),qT(t)). 

Under  Assumptions  A.l  - A. 4,  we  will  construct  an  asymptotic  solution  which 
is  valid  in  a region  away  from  the  end  points.  We  will  substitute  the  series  (3.5) 
and-  (3.8)  into  the  full  problem  (2.23)  - (2.27).  The  resulting  equations  may  be 
differentiated  j times  with  respect  to  e.  By  setting  e = 0,  we  obtain 


For  j = 0, 


xo  = f{xZ,z%,P%,q£)  (3.9) 

Po  = d{xo  izo  iPo  i Vo)  (3.10) 

0 = F(x%,z°,p%,q£)  (3.H) 

0 = G(x%,zg,p%,q%)  (3.12) 

0 = cf>t  + H(x°,z%,p%,q°)\t=To.  (3.13) 

For  j > 1, 

x°  = f*x°  + fPP°  + f.zf  + fqq°  + rg  (3.14) 

P?  = 9*x°  + 9pP°  + 9zZ°  + gqqf  + rg 


(3.15) 
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0 = Fxx°  + FpP°  + F2z°  + Fqq°  + r% 

(3.1G) 

0 = Gxx°  + Gpp°  -\-Gxz°  + G,t‘\i  + 

(3.17) 

0 = r°Ti  + H,x°  + Hrp°  + T°  |1=T> 

(3.18) 

where  r°  for  i = 1 • • • 6 are  successively  known  functions. 

Now,  we  will  derive  the  boundary  conditions  for  the  problem  (3.9)  - (3.13)  and 
(3.14)  - (3.18).  Assuming  that  xR(a,  e)  is  negligible  at  t = 0,  substitution  of  (3.1) 
into  the  initial  condition  of  the  slow  state  of  the  full  problem  gives 

OO 

I>?(0)eJ  + xJL(0)eJ+1)  = x0. 
j=o 

Matching  the  coefficients  of  equal  powers  of  c gives 


x°(0)  = x0 

(3.19) 

x°(0)  = -x^,(0)  j > 1. 

(3.20) 

Also  assuming  that  xL(r,e)  is  negligible  at  t = T,  substitution  of  (3.1)  and 
(3.8)  into  the  terminal  condition  of  the  slow  state  of  the  full  problem  gives 

OO 

^2(xf(To  + ePi  + e2T2  + • • •)(*  + x^(0)eJ+1)  = Xj- 
j= o 

By  matching  the  coefficients  of  equal  powers  of  e 

(To)  = %T 

(3.21) 

x°(T„)  = -xf_1(0)  + r°r,  +r°  j>l 

(3.22) 

where  r°-  and  r are  successively  known  values. 

Definition  3.1.  The  asymptotic  series  solution  x°(t),p°(t),  z°(t), 
is  called  the  jth  order  outer  solution  of  the  full  problem  if  it  satisfies 

qf(t)  and  Tj 
the  recursive 

2 n dimensional  TPBVP  (3.9)  - (3.18)  with  the  boundary  conditions  (3.19)  - (3.22). 
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We  remark  that  the  zeroth  order  outer  solution  {x°  ,p°  ,z%  ,q°  ,T0)  is  same  as  the 
solution  of  the  reduced  problem.  Once  the  zeroth  order  outer  solution  is  obtained, 
the  higher  order  solution  can  be  obtained  by  solving  the  linear  inhomogeneous  time 
varying  TPBVP.  Since  B(t)  is  invertible,  zf(t ) and  q° (t)  can  be  uniquely  deter- 
mined in  terms  of  xf(t)  and  p°(t)  from  the  algebraic  equation  (3.16)  and  (3.17). 
Assumptions  A. 2 and  A. 3 will  guarantee  the  existence  and  uniqueness  of  the  lin- 
ear inhomogeneous  time  varying  TPBVP.  One  powerful  method  to  solve  the  linear 
inhomogeneous  TPBVP  is  the  backward  sweep  method  [25]. 

3.3  Boundary  Layer  Expansion 

The  outer  solution  constructed  in  the  previous  section  does  not  satisfy  all  the 
boundary  conditions.  In  order  to  obtain  a uniformly  valid  solution  throughout  the 
entire  interval,  boundary  layer  correction  terms  known  as  the  left  boundary  layer 
solution  and  the  right  boundary  layer  solution  will  be  necessary.  Assuming  that  the 
right  boundary  layer  terms  are  negligible  in  the  left  boundary  layer,  from  (3.1)  - 
(3.4)  the  left  boundary  layer  terms  must  satisfy  in  time  scale  r 

dxL(r,  e) 
dr 

dpL(r,e) 
dr 

dzL(r,e) 

dr 

dqL(r,e) 
dr 

with  the  initial  condition 


'dx(t,e)  dx°(t,e) 


= e 


= e 


dt 

dt  j 

dp(t,e) 

dp°(t,c)\ 

dt 

dt  J 

( dz(t,e) 

dz°{t , e) 

V dt 

dt 

( dq{t,t) 

dq°  (t,  e) 

\ dt 


dt 


t — CT 


t—CT 


t — €T 


t — tT 


(3.23) 

(3.24) 

(3.25) 

(3.26) 


2l(0,  e)  4-  z°{ 0,  t)  = 20. 


(3.27) 
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We  may  differentiate  (3.23)  - (3.27)  j times  with  respect  to  e and  set  t — 0. 
Then  the  resulting  equations  will  be  given  by 


For  j = 0, 


dr 


= /(4\4  + 4(4, 4,4  + 4(4) 


dpo(r) 

dr 


<*4(4 

dr 

<*4(4 

dr 

4(0) 


-/(4\4,4,4) 

$(4,4  + 4(4, 4, 4 + 4(4) 
-5(4,4,P?,4) 

44,4  + 4(4, 4, 4 + 4(4) 
G(4  , 4 + 4(4, 4, 4 +4(4) 
*0  - 4(0) 


For  j > 1, 

<*4(4 

dr 

<*4(4 

dr 

■PM. 

dr 

<*4(4 

dr 


fz{r)zf(r)  + /,(r)4(r)  + rf^r) 

gz(r)zf(r)  + £,(44(4  + 4j(4 

F*(T)Zj(r)  + F„(T)qf(T)  + 4(r) 

G*(r)4  (T)  + <4(44(4  + 4,(4 


where  for  example, 

/«(t)  = ^(x„°(0),zo°(0)  + 4(r),P?(0),9„°(0)  + ?0i(T)) 


(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 

(3.36) 

(3.37) 


7+(r)  are  successively  known  functions  for  i = 1 • • -4. 
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Definition  3.2.  The  asymptotic  series  solution  Xj(t),  and  qf(r)  is 

called  the  jth  order  left  boundary  layer  solution  of  the  full  problem  if  it  satisfies 
the  recursive  differential  equations  (3.28)  - (3.37).  Moreover,  we  require  that  there 
exist  some  positive  constants  C and  v such  that  for  all  r 6 [0,  oo), 

K%)l  + \Pj(r)\  + \zf(r)\  + \qf(r)\  < Ce~” . 


We  remark  that  only  the  initial  conditions  of  the  fast  state  (3.32)  and  (3.37) 
are  imposed  on  the  differential  equations.  In  order  to  solve  the  left  boundary  layer 
problem,  q^( 0)  for  j = 0, 1,2,  • • • must  be  so  chosen  that  the  resulting  solution  zf(r) 
and  qj(r)  satisfying  (3.30)  - (3.31)  or  (3.35)  - (3.36)  decay  exponentially  to  zero 
as  t goes  to  infinity.  Then,  xj'(r)  and  pf(r)  can  be  obtained  by  simple  integration 

. dx^(r)  , dpz(r)  ... 

smce  — and  are  known  functions  of  a stretched  time  scale  r.  In  order  to 

guarantee  the  existence  of  the  left  boundary  layer  solution  of  the  full  problem,  we 
need  the  following  assumption. 

Assumption. 

A. 5-  There  exists  a 2 m x 2 m continuously  differentiable  matrix  M(t)  for  all 

0 < t < Tq  such  that 
1. 


2. 


B\{t)  0 

0 


A/„(0)  # 0 A/22(r„)  j'  0 


3.  Bi(t)  has  m eigenvalues  with  negative  real  parts  for  each  t,  0 < t < T0. 
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Under  Assumptions  A.l  and  A. 5,  it  was  proved  by  Freedman  [22]  that  for 
a sufficiently  small  positive  constant  6 , the  left  boundary  layer  equation  (3.28)  - 
(3.37)  has  a unique  solution  decaying  exponentially  to  zero  as  r goes  to  infinity 
for  each  j = 0, 1 , 2 • • * if  |z^(0)|  < 6.  Here,  a smallness  condition  of  | ( 0 ) | is 
inevitable  in  a general  nonlinear  context.  However,  this  condition  may  be  eliminated 
in  some  cases,  for  example,  quasi-linear  systems.  The  procedure  to  construct  the 
left  boundary  layer  solution  stated  in  this  section  is  the  same  as  in  a fixed  final  time 
case.  Even  if  the  existence  of  the  left  boundary  layer  solution  is  guaranteed,  it  may 
be  prohibitively  difficult  to  solve  the  left  boundary  layer  equation. 

We  will  proceed  to  construct  the  right  boundary  layer  solution.  Similarly  as 
in  the  left  boundary  layer,  in  a time  scale  the  right  boundary  layer  solution  must 
satisfy 


dxR(a, 

<0 

da 

dpR(a, 

<0 

da 

dzR(a , 

<0 

da 

dqR(a , 

0 

da 


( dx(t,e)  dx°(t,e) 

\ dt  dt 


t — To+cTi  +e^T2-f** — to 


( Mt, e) 

V dt  + 


t=To  +eXi  +c2T2-\ e<r 


e 


( dz(t,e ) 

V dt 


+ 


dz°(t,e)\ 
dt  J 


t=To  +(.T\  +e2  Ti  H f a 


<*g(*>c)  , dq°(t,eY 
dt  dt 


t—To  +(.T\ +e2  T2  H e<r 


with  the  initial  condition 


(3.38) 

(3.39) 

(3.40) 

(3.41) 


9fi(0,  e)  + q°  (T„  + t r,  + e2T2  + ••-,£)  = 0.  (3.42) 

We  may  differentiate  (3.38)  - (3.42)  j times  with  respect  to  ( and  set  f --  0 . 
Then  the  resulting  equations  will  be  given  by 
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For  j = 0, 


For  j > 1, 


u 


dx^{a) 

da 

= -/(*?, *0°  + <iS  + rfW)) 

(3.43) 

dpo(<r) 

da 

= + *o(<0>Po>9o  +9o(<7)) 

5 jPo  5 % ) 

(3.44) 

dzg(a) 

da 

= -F(xo,zo  + 2u(<T)>Po>9o  +g£V)) 

(3.45) 

dq*(a) 

da 

= -<?(*?,  *o  + 2*(oO>Po>9o  +9^(cr)) 

+G(x°,Zo°,p°,?f) 

(3.46) 

«?(  o) 

= -?o°(0) 

(3.47) 

dxf(a) 

da 

= -/.(*)*?(*)  - /»??  (*)  + »i» 

(3.48) 

dpf{a ) 
da 

= -9»{?)zf  M - 9q{°)(l?(v)  + rj(<r) 

(3.49) 

da 

= -F.(<,)zf(a)  - F,W)qf(c)  + r*(<x) 

(3.50) 

dqf{a) 

da 

= -G»zf  (<r)  - G,(a)gf  (a)  + r4» 

(3.51) 

<7?(0) 

= -9f(r0)  + rg 

(3.52) 

xessively  known  functions  for  i = 1 ■ • • 5 and 

df( 

°(r0),^(T0)  + z«(a),p°(T0),q°(T0)  + q*(a)). 
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Definition  3.3.  The  asymptotic  series  solution  xf(o),  p^(a),zf(a)  and  qf(a)  is 
called  the  jth  order  right  boundary  layer  solution  of  the  full  problem  if  it  satisfies 
the  recursive  differential  equations  (3.43)  — (3.52).  Moreover,  we  require  that  there 
exist  some  positive  constants  C and  v such  that  for  all  a 6 [0,  oo), 

l*jV)l  + |p?HI  + \zf(o)\  + |9f  (a)\  < Ce~ua . 

As  in  the  left  boundary  layer,  only  the  initial  conditions  of  the  fast  costate 
(3.47)  and  (3.52)  are  imposed  on  the  differential  equations.  In  order  to  solve  the 
right  boundary  layer  problem,  zp( 0)  for  j = 0, 1, 2,  • • • must  be  chosen  such  that  the 
resulting  solutions  and  <7^(<t)  decay  exponentially  to  zero  as  <r  goes  to  infinity. 

Under  Assumptions  A.l  and  A. 5,  the  right  boundary  layer  problem  has  a unique 
right  boundary  layer  solution. 

3.4  Discussions 

In  summary,  the  procedure  to  obtain  an  asymptotic  power  senes  solution  is 
given  by 

1.  From  the  reduced  problem  (2.33)  - (2.38),  we  obtain  a unique  reduced  solution 
xr{t)iPr{t)i  zr(t),  qr(t)  and  Tr. 

2.  Based  on  the  reduced  solution,  we  obtain  z£(r),  q% (t),z§ (a)  and  qp(a)  from 
(3.30)  - (3.32)  and  (3.45)  - (3.47). 

3.  Since  the  zeroth  left  and  right  boundary  layer  solution  are  known,  Xq(0)  and 

(0)  are  computed  from  (3.28)  and  (3.43). 

*4(0)  = f°^dr 

J oo  U.T 

*4(0)  = f^da 

J oo  d<7 
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4.  Since  we  can  determine  x?(0)  and  x?(T0)  from  (3.20)  and  (3.22),  the  first 
order  outer  solution  can  be  obtained  from  (3.14)  - (3.18). 

5.  Proceeding  recursively,  the  j th  order  outer  solution  allows  us  to  solve  the  j'th 
order  left  and  right  boundary  layer  solution  and  the  boundary  condition  of  the 
j + 1 st  order  outer  problem  can  be  obtained  from  the  j th  order  left  and  right 
boundary  layer  solution. 

After  determining  the  power  series  solution  in  each  time  scale  up  to  the  nth 
order,  an  approximate  solution  of  the  full  problem  can  be  made  as  a sum  of  the 
power  series  solution  from  (3.1)  - (3.4). 


CHAPTER  IV 
ASYMPTOTIC  VALIDITY 


In  the  previous  chapter,  we  constructed  an  asymptotic  series  solution  in  order 


to  obtain  an  approximate  solution  of  the  full  problem.  In  this  chapter,  we  will 
show  that  the  construction  method  described  in  the  previous  chapter  will  give  a 
good  approximate  solution  of  the  full  problem  if  e and  the  size  of  the  zeroth  order 
boundary  layer  jumps  in  the  fast  variables  are  sufficiently  small. 

It  is  convenient  to  make  the  following  abbreviations: 


2/(0  = 

x(t) 

1 

C-+. 

» 

6^ 

II 

. ?(o . 

<7(0  . 

Thus,  y(t)  represents  the  slow  state  and  the  slow  costate  and  u;(7)  represents  the 
fast  state  and  the  fast  costate. 


For  any  integer  N > 0,  we  define  the  partial  sums  by 


Vn 

£f=O2/?(0^  ' 

. . 

1 

■«-» 

MJ 

°sT 

4 

w 

■ 

Vn 

IS  Iff  (r)e»‘ 

. . 

. E"o  «f(r)e»‘  _ 

Vn 

EjLo1  yf  (*)* 

. . 

. EjlotofOry 

Vn 

2/w(0  + e2 In(t)  + ey^(o-) 

wN(t)  + Wft(r)  + Wft(cr) 

(4.1) 


(4.2) 


(4.3) 


(4.4) 
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Tn  = Er/  (4.5) 

i= o 

Let  y(t,e),  w(t,e)  and  T(e)  be  the  solution  of  the  full  problem.  Then  y(t,e),  iv(t,e) 
and  T’(e)  can  be  expressed  by 


y(t,e) 

VN(t ) 

y'(t) 

= 

+ eN 

_ U>(M)  _ 

wN(t) 

T(e)  = TN  + eNY 

T(e)  -t  _TN  + eNY  - t 

t t 

where  eNy',  tN  w*  and  cNj*  denote  the  remainder  terms. 

The  objective  in  this  chapter  is  to  prove  that  i)  the  solution  of  the  full  problem 
y(M),  w(t,  e)  and  T(e)  exists  under  Assumptions  A.l  - A. 5 described  in  Chapter  II 
and  ii)  y(t,e),  w(t,  e)  and  T(t)  approach  the  asymptotic  solution  yN,  wN  and  T/v  as 
eN+1  goes  to  zero.  The  principal  tool  to  be  used  in  the  proof  is  a Banach  space  fixed 
point  theorem  [26]. 

4.1  Nonlinear  TPBVP  for  the  Remainder  Terms 

We  will  derive  a nonlinear  TPBVP  governing  the  remainder  terms  ym , wm  and  7*. 
One  can  see  that  for  sufficiently  small  e there  exists  T > Tn  such  that  Assumption 
A. 5 also  holds  for  an  interval  [0,T],  Thus,  it  will  be  assumed  that  Bx(t ) defined  in 
Assumption  A.5  has  m eigenvalues  with  real  parts  < -2 v < 0 for  each  t e [0,T], 
We  also  assume  that  the  outer  solution  is  extendable  such  that  the  outer  solution  is 
defined  throughout  the  compact  interval  [0,T]. 
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Using  (4.4)  and  differentiation  of  (4.6)  shows  that  the  remainder  terms  must 
satisfy 


iP  , dVs  dy^  N 

+ e y - u(y,  w) 


Vn  + 


dr  da 


■-o  . dwLN  dw%  N+l 

h e ^ w = U(y,w) 


ewN  + 


dr  da 


(4.7) 

(4.8) 


where 


u(y,w)  = 


dt 


U{y,w ) = e 


dw(t , t) 
dt 


f{x,z,p,  q ) 

g(x,z,p,  q) 

F(x,z,p,q) 

G(x,z,p,q) 


By  the  construction  of  the  outer  solution  in  the  outer  problem  (3.9)  - (3.18), 
y )y  and  w1^  must  satisfy 


■-O 

Un 


= (u(yZ,™Z)) 


N 


= «($,«#) + 0(e"+1) 

where  (• )N  represents  all  the  terms  containing  ek  for  0 < k < N. 
Similarly, 


(4.9) 


>-0  -o 


N+ 1> 


(4.10) 


We  introduce  the  size  of  the  jth  order  boundary  layer  jumps  8j  such  that 


Also,  we  define  8 by 


”4(0)1,  |<(0)|<«,  0 <;<«. 


8 = max  8j. 

0<j<N  J 
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From  the  construction  of  the  left  boundary  layer  solution  in  the  left  boundary 
layer  problem  (3.28)  - (3.37), 

“|r  = + + 

= ^ + ^,^d-u>^)-u(^,u}°)  + 0(e^e-"T)  (4.11) 

^ = {U(C  + eyi,wZ  + ^i)~U(yZ,wZ))N 

= U{y°N^ty^,w^  + wLN)-U{y^,w^)-VO(tN+l8e-UT).  (4.12) 

In  the  right  boundary  layer,  the  construction  (3.38)  - (3.41)  gives 

dyR 

= - ecr)  + £Vni  ™n{Tn  — ecr)  + wR) 

+u(y°(TN  - ecr),  w%{Tn  - e<T)))^v_1 
= (~u(yN(t  ~ tNl*)  + ey%,  w%(t  - eN~f*)  + u$) 
+u(y°(t-eN7*),w°(t- 
= ~u(yN (t  — eN7*)  + q/$)  Wn(t  — eN7*)  + wR) 

+u(y%(t  ~ eV),  w°N(t  - eN 7*))  + 0(eN8e~^)  (4.13) 

dwRr 

~&T  = + + 

+£/($(*  - eV),4(f  - eV))  + 0{tN+x8e-v°).  (4.14) 

Substitution  of  (4.9)  - (4.14)  into  (4.7)  and  (4.8)  gives 

eNym  = u(y%  + ey ^ + eyR  + eNy*,w%  + + wR  + eNw‘) 

u(y n + ey^,  w%  + wfr)  + u(y%(t  — eN7*),  w%(t  — eN7*)) 

~u(y%(t  - eV)  + ty%,  W%(t  - e V)  + u>£) 
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eN+'w* 


+0(tN+1 ) + 0(tN  8e~l'T)  + 0(eNSe-w) 

U{Vn  + + tyR  + eNym,u>H  + + Wm  + eNw *) 


(4.15) 


U(Vn  + eVN^wN  + wn ) + U(y%(t  — eN  1 *),w°(t  — eN  7*)) 

-U{y%(t  - eN7*)  + «7*,u>£(i  - £N  1*)  + ™n) 

+0{eN+l)  + 0(eN+18e~l/T)  + 0(eN+18e-'cr).  (4.16) 


We  rewrite  the  above  equations  by  splitting  the  right  hand  sides  of  (4.15)  and 


(4.16)  into  linear  parts  and  the  remaining  nonlinear  parts. 


V * = uyV * + uww*  + u*(y*,w *,  7*,  e)  (4.17) 

ew*  = Uyy*  + Uww*  + U*(y*,w *,  7*,  e)  (4.18) 


where 


,7V  * / * * * \ 

e u (y  ,w  , 7 ,c) 


u(yjv  + ey^  + ey$  + ^ + eNw*) 

~u(Vn  + w%  + wln)  + u(y%(t  - e V)>  w°(t  - £N7*)) 
~u(y%(t  ~ eN-y*)  + ey%,  w°N(t  - ew7*)  + w%) 

—tNuyy*  - eNuww* 


+0(eN+1)  + 0(eN8e-'T)  + 0{tN8e~ua)  (4.19) 

iu*,7*,e)  = U{y°N  + ey^  + ty^  + tNy\w^  + w^^w^  + tNw') 

~^{yN  + ey x,  + wfr)  + U (y°(t  — £^7 *),  u>%(t  — £^7*)) 

~U(y°(t  - tN-y*)  + «/$,«>£(*  - £^7*)  + wjv) 


-^y*  - 
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+0(e/v+1)  + 0(eN+l6e~l/T)  + 0(eN+lSe~ u°) 


(4.20) 


where  uy  is  defined  by 


9u ,Q  -o  \ 
uv  = Qjjil /*>«#»«) 


and  Uy  and  Uw  are  also  similarly  defined. 

The  boundary  conditions  for  the  remainder  terms  can  be  determined  as  follows. 
The  initial  x is  given  by 


*o  = xjv(O)  + eNx*(0) 


= xjv(O)  + ex^(O)  + «£( 
Since  x^(0)  + cx^(O)  = x0  by  construction, 


-R/Tn  + t 7* 


) + cNx*(0). 


i*(o)  = 


(4.21) 


The  boundary  condition  at  the  terminal  time  is  obtained  by  considering  the 
final  x. 


XT  — x^(T/v  + eyV7*)  + exN( 


N^.*\  , -l /Tn  + eN 7* 


) + ex*(0)  + eNx*(TN  + eV) 


xn(Tn)  -fi  ex^(O)  + ex^( 


_L.TN  + eNY 


) 


+ t1"*"*  , ,N 

JT„ 


x„{s)ds  + eNx*(TN  + eN7*). 


Since  xn(Tn)  + ex^(O)  = xj  + 0(eN+1)  by  construction, 


*t rp  , n *\  -N+i  -L  / Tn  T e^7*  \ 1 [Tn+*nY  ,q 

x (Tn  + e 7 ) _ — e xN(  ) — xN(s)ds  + 0(e).  (4.22) 

£ 6 JTn 


Similarly, 


*•(  0)  = -<-"*»(  Zk±iV) 


(4.23) 


*,rp  AT  *\  -N-LfTtf  + e1*  7*  1 rTN+cN-y- 

<7  (Tn  f£  7 ) = -e  qN( ) ^ /_  9&(s)ds  + 0(e)(4.24) 

e e ./tv 
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The  transversality  condition  becomes 


0 = + H(yN(t)  + t‘vym(t),wN(t)  + eNw*(t)) 


t=T\r+eN7* 


(4.25) 


We  also  rewrite  (4.25)  into  linear  parts  and  the  remaining  nonlinear  parts. 


Cfyi(7w)  + C2  y2{TN)  + C3j‘  = -irty\iv-,;\  0 


(4.26) 


where 


Cx 

C2 

eNH‘(y\w\r}e) 


Hyi(y°(TN),w°(TN)) 

Hy,(y%(TN),w%(TN)) 

MTn ) + C*y\°(TN)  + CU2°N{TN) 


<M0  + H(yN(t)  + eNy*(t),wN(t)  + eNw*(t)) 


t=TN+iN-y* 


-tNC*y\(T„)  - tNCly-2(TN)  - eNC3l\ 


(4.27) 


In  order  to  obtain  some  necessary  properties  for  H*(y*,  w*,  7*,  e),  we  need  the 
following  lemma. 

Lemma  4.1.  For  the  partial  sums  (4.4)  and  (4.5)  constructed  in  Chapter  III, 

M^n)  + H(y%(TN)  + ey^(~)  + ey%(0),w%(TN)  + u>&(^)  + tt$(0)) 

= 0(eN+1). 

Proof-  Since  y^ijf),  are  0(eN+1), 

MTn ) + H(y^(TN)  + ty^{-T-)  + ey%(0),w%(TN)  + ^w(~)  + u>$(0)) 

= MTn)  + H(y°(TN)  + ey*(0),  w°(TN)  + u$(0))  + O^1).  (4.28) 
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Now  consider  a function  V(a)  defined  by 

V(cr)  = H(y%(TN  - ea)  + ey*(a),w%(TN  - ea)  + wfi(a)) 

~H {Vn{Tn  - ecr),  w%(TN  - ea)). 

Since  y*(^f ) and  w^(^f)  are  0(e;v+1), 

V(~)  = h(v°n(  0) + ^(v)’”«(°) + "»(—)) 

= 0(eiV+1). 


Using  (4.9)  (4.10)  and  (4.13)  (4.14),  differentiation  of  V(a)  with  respect  to  a 


gives 


Hence,  for  some  K , 


By  integrating, 


dV(a) 

da 


= 0(eN+1Se~l'a). 


dV  (a) 


da 


< KeN+16e~l/<T . 


< 

VftlL) 

+ 

e 

jZu-  da 
€ 

= 0{eN+1). 

From  the  construction  of  the  outer  expansion  (3.18), 

MTn)  + H(y°(fN),w°(TN))  = 0(eN+1). 

Hence, 


<t>t{TN)  + H(y%(TN)  + ey^(0),w^(Tiv)  + u^(0)) 
= MTn)  + H(y°(TN),w°(TN))  + U(0) 

= 0{eN+1). 


(4.29) 
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By  combining  (4.29)  with  (4.28),  the  proof  is  complete.  ■ 

For  the  nonlinear  parts  of  (4.17)  and  (4.18)  which  will  be  considered  as  forcing 
terms,  we  state  some  properties  to  be  useful  in  proving  our  main  result. 

Lemma  4.2.  Let  = e"^e  + c-^n-«jv+1  A'-OA  Given  I<  > 0,  u\y\ w*,  7*,  e), 
U*(y',w\  7V)  and  H*(y*  , re*,7*,c)  have  the  following  properties. 

1.  There  exist  a kx  > 0 and  e1  > 0 such  that  when  |y*|,  |iu*|,  |7*|  < AT  for 

0 < e < ei 

!)  K(yW,7*,e)|  < kx  (e  + 8xl>  + (e  + 80^)(\y*\  + |u>*|)  + (e  + 60)^|7*|j 

2)  \U*(y*,  w*,  7*,  e)|  < kx  (e  + (e  + 60)(|y*|  + \w*\  + |7*|)) 

. 3)  I H (y  ,w  ,7*,e)|  < kx  (e  + (e  + ^o)(|l/*|  + |tw*|  + |7*|))  • 

2.  For  each  k2  > 0,  there  exists  a A^fcj)  such  that  |y|,  |y|,  |u>|,  |i2>|,  |-y|,  |7|  < AT, 
0<£o<  Ai(A:2)  and  0 < e < Xx(k2)  imply 

1)  |w*(y,ih,7,e)  - u*(y,w,  7,e)|  < k2  - y\  + |u>  - w\  + - 7|j 

2)  |t/*(y,u),7,e)  - U*(y,  w,  7,  e)|  < k2  (|y  - y|  + |tu  - w\  + |7  - 7I) 

3)  \H*(y,w,j,  e)  - H*(y,u>,  7,  e)|  < k2(\y  - y\  + |ru  - tu|  -f  |7  - 7I)  . 

ProoL  Using  0{e~utll)  = 0(ek)  for  t > TN/2  and  0(e-^Tlt~t^e)  = 0(ek)  for 
t < Tn/2  for  an  arbitrary  positive  k,  one  can  show  that  in  a ball  |i/*|,  |re*|,  |7*|  < AT, 

«*(0, 0,  0,  e)  = 0(e  + 6t/>)  U*( 0,  0, 0,  e)  = 0(e) 

uy*(y  ,w  >7*>  e)  = 0(e  + Soip)  u* . (2/*,  iw*,7*,  e)  = 0(e  + 60^) 


= O(e  + £oV0  U*.(y*,w',  7*,e)  = 0(e  + <W>) 
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ui  m(y  ’ w >7*,e)  - 0(e  + 60  )0(^)  U*.(y*,  wm , 7*,  e)  = 0(ct/>  + 6O0). 

Hence,  properties  1.1),  1.2),  2.1)  and  2.2)  directly  follow  from  the  mean  value  the- 
orem. 

Since  = O{e  + 60)  from  (4.27), 

\H*(y,w,j,e)  to, 7,  e) | < fc(e  + <50)|7-7|  (4.30) 

for  some  A;  > 0. 

For  simplicity,  we  define 

7//(0  = /7yT  (MO  + = HZ(y°(TN),w°(TN)) 

H l = Hy(VN{t)  + eNy(t),wN(t)  + eNw(t)), 

and  H*,,  H * and  II sw  are  also  similarly  defined. 

In  order  to  estimate  \H*(y,  w,  7,  e)  — H*(y,  w,  7,  e)|,  we  rewrite  (4.27)  into  the 
following  form, 

ir  = hi  + h; 

where 

Hi  = ^(Vn{Tn)  + eNy*(TN),wN(TN)  + eNw*(TN))  - eNH°y*(TN) 

+MTN  + eV)  - tNMTN)r 

[Tn+cn 7*  , 

H 1 = JTn  (- Hl{VN  + £Nr)  + Hl{iuN  + eNw*)-H°y°(TN))ds . 

Since 

Hly.  = H'(TN)-H°  = O(e  + 60) 


H*lw.  = Hl(TN) 
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— Hwiy^iTw),  w^(T^))  + 0(e  + So) 

= O(e  + 60 ) from  (4.10), 
one  can  deduce  that  for  some  k > 0, 

\Hl(y,w,i,e)  - Hi(y,w,i,e)\  < k(e  + S0)(\y  - y\  + |u>  - u>|). 
\H;(y,w,i,e)  — H^y,™,  7,e)|  can  be  expressed  by 

\HZ(y,w,l,z)  ~ HZ(y,w,i,e)\ 

< L I - H’\\tN  + e"y|  + eN\ H'\\y  - jj\ds 
J In 

rTjv+e^  7 A „ . v 

+ L I HL  - H1\\™n  + tNw I + tN\Hsw\\w  - w\ds. 

JTn 

We  use  the  fact  that  for  some  k > 0, 

=0(1)  Hy  ~ O(l) 

ti,v  + eN<i>  = 0(i±M)  ///  = 0(£  + W) 

and 

l^--ffy|  <fc|y-y|  |y-y|<fc|y-y| 

|//£  — //^|  < A:|u)  — uij  |u>  — rh|  < — |u)  — th|. 

Hence,  equation  (4.32)  becomes 

- H;(y,w,j,t)\  < k(e  + <5o)(|y  — y|  + |ti>  - ii|) 


(4.31) 


(4.32) 


(4.33) 


for  some  k > 0. 
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By  combining  (4.30)  and  (4.31)  with  (4.33),  property  2.3)  is  established.  By 
using  (4.27)  and  Lemma  4.1, 

i/*(0,0,0,e)  = (^(TV)  + II(yN,  tu/v))  = 0(e). 

Hence, 

\H*(y* ,w* ,'y*,  e)|  < |77*(0, 0,  0,  e)|  + u>*,  7*,  e)  — //*(0, 0,  0,  e)| 

< fci(e  + (c  + 60)(|y*|  + |w*|  + l7*l))- 

This  completes  the  proof.  ■ 

4.2  Transformation  of  the  TPBVP 

Now,  with  the  aid  of  Lemma  4.3  taken  from  P.  Ilabets’s  work  [27],  we  will 
transform  the  system  (4.17)  and  (4.18)  into  a tractable  form  such  that  the  slow 
variables  and  the  fast  variables  are  weakly  coupled. 

Lemma  4.3.  Suppose  that  i)  Fi(t),  F2(t),  Gi(t)  and  G2(t)  are  continuous  and 
ii)  G2{t)  has  m eigenvalues  with  real  parts  < —2u  < 0 and  m eigenvalues  with  real 
part  > 2u  > 0 for  each  t 6 [0,T],  where  Fi(t),  F2(t),  Gi(t)  and  G2(t)  are  2 n x 2 n, 
2 n x 2m,  2m  x 2 n and  2m  x 2m  dimensional  matrices,  respectively.  There  exists 
e2  > 0 such  that  for  any  e G [0,  e2]  the  equations 

eR{t)  = G2(t)R(t)  - eR^F^t)  + eR(t)F2(t)R(t)  - Gtit) 

eS(t)  = e(Fl(t)-F2(t)R(t))S(t)-S(t)(G2(t)  + eR(t)F2(t))-F2(t) 
have  solutions 

R=  R(t,e)=  G2  1Gi  + 0(e  + ip)  where  xj)  = e~ut^  + 


S = S(t,e)=  0(1). 
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Let 


F\  — Uy  i*2  — UW  G\  — Uy  G2  — Uw 

Since  the  conditions  for  Lemma  4.3  are  satisfied,  R(t ) and  S(t)  have  solutions 

R(t)  = U~1Uy  + 0(c  + rP)  S(t)  = 0{  1). 

Consider  the  change  of  variables, 


a 


P 


I + eSR  tS 
M~lR  A/'1 


“I 

* 

y * 

1 

* 

3 

1 

(4.34) 


(4.35) 


where  M{t)  is  defined  in  Assumption  A. 5. 

Lemma  4.4.  With  the  change  of  variables  (4.34),  the  TPBVP  (4.17)  and  (4.18) 
can  be  transformed  into 


y* 

I -eSM 

a 

w* 

-R  M + eSRM 

J _ 

a = A(t)a  + /*(a,/?,7*,e) 

(4.36) 

eP  = B(t)f3  + m*(a,p,  7*,e) 

(4.37) 

where 


A(t)  = Uy-UwU-lUy 


B(t)  = M~XUWM 


Bi(t)  0 

0 -BT(t) 


Moreover,  l*(a,(3, 7*,e),  m*(a,  /3, 7*,  e)  and  7*,e)  have  the  same  properties 

as  in  Lemma  4.2.  Thus, 
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1.  There  exist  a kx  > 0 and  e,  > 0 such  that  when  |a|,  \/3\,  |7*|  < Kt  for 

0 < e < 

1)  K (a5  /?>  7 ? e)l  5:  + <50  + (e  + + \P\)  + (e  + S0)  — 17*|^ 

2)  | m (a,  /3, 7*,  e)|  < ki  (e  + ( e + «50 ) ( | o; | + \/3\  + |7* | )) 

3)  \H-(a,/3,r,  e)\  < *a  (e  + (e  + 50)(|a|  + |/5|  + |7*|)) . 

2.  For  each  k2  > 0,  there  exists  a A i(k2)  such  that  |a|,  |a|,  |/?|,  \j3\,  |-y |,  |^y|  < K e, 
0 < So  < Xi(k2)  and  0 < e < Ai(&2)  imply 

1)  |!*(d,M,e)  - C(a,/9,7,£)|  < k2  ^|d  - d|  + |/j  - 0\  + ^7  - A 

2)  |m*(a,^,7,c)-m*(a,^,7,e)|<fc2(|a  — a|  + |^  — ^1  + 17  — 7|) 

3)  |//*(a,^,7,e)  - H*(a,  ^,7,  e)|  < k2  (|a  - d|  + \0  - j}\  + I7  - 7I)  . 

Proof.  By  substituting  (4.34)  into  (4.17)  and  (4.18),  one  can  obtain  (4.36)  and 

(4.37).  Also,  l*(a,  /?, 7*,  e)  and  m*(a,  (3, 7*,  e)  are  expressed  by 

= Aa+Af3  + u* 
m*(a,/?,7*,e)  = B/3  + LT 

where 

A = uw  (1 U-'Uy  -R)  = 0(e  + xj)) 

A = e2  (uySSRM  — uwRS SRM  — uySRSM  + uwRS RS Al) 


= 0(e2) 
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um 

I + eSR 

S 

u * 

U* 

tM~'R 

M-1 

u- 

D = e[M  lRuwM  + M~XUWSRM  - M~x MM~lUwRSM') 
+e2  (m~xRuwSRM  + M~lMM~lRSM 


- M 1MM~1SRM  - M~'RuwRSM ) = 0(e). 


From  (4.34),  for  some  K > 0, 

M < |y*|  + eK(\y*  \ + |u>*|) 

\P\  < *(|y*l  + M)- 


Hence,  the  properties  of  /*(a, /?,  7*,  e),  m*(a, /?,  7*,  e)  and  //*(a,  /?, 7*,  e)  follow  di- 
rectly from  the  definition  of  /*(a,  /?,  e),  m*(a,  /?,  e)  and  //*(a,  /?,  7*,  e).  ■ 

We  remark  that  H*(a,  (3, 7*,  e)  in  Lemma  4.4  is  obtained  by  substituting  (4.34) 
into  e). 

In  terms  of  the  new  variables,  the  boundary  conditions  and  the  transversality 
condition  will  be  derived.  Let 


- 

- 

r 

^nXn  0 

P2  = 

0 0 

Q = 

X m 0 

1 

O 

O 

IjiXn  0 

1 

O 

O 

By  using  (4.35),  the  boundary  conditions  (4.21)  - (4.24)  and  the  transversality 
condition  (4.26)  become 


ai(0) 

0 

_ ai(fjV  + eNY) 

1 rfAr+eN7*  -Or  , 

tN  JTn  xN(s)ds 

+ 0) 

+eP2SMp(fN  + eNY)  + 0(e-‘/^+£^*)/£) 


(4.38) 
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M(7V) 


A(0) 

HTn  + eNr) 


= QRa( 0)  + (/  - Q)Ra(TN  + e^*) 


0 

1 rTW+t'V  =-0 


r 

JIT  if. 


qN(s)ds 


+ 0(e~,'(TN+tNl">/c) 


(4.39) 


Ci<Xi(Tn)  + Cja2(T)v ) + C3 7*  = 7*,  e) 


+e  [cf  Cf]  (A5M  + P2SA/)  /?(7W) 


(4.40) 


where 


M(7*,  c)  = Q(I  + eS.R)M (0)  + (/  - Q)(I  + tSR)M{TN  + eNf).  (4.41) 

Consider  the  linear  parts  of  (4.36)  and  (4.37).  Let  ©i(f,s)  and  02(t,s)  be  the 
fundamental  matrices  satisfying  the  following  matrix  differential  equations 

0i(C-s)  = s)  01(s,s)  = /mxm 

02(M)  = --Bf(t)Q2(t,s)  02(s,s)  = ImXm. 

It  was  shown  by  Flatto  and  Levinson  [28]  that  0i(<,s)  and  02(*,s)  have  the 
following  bound  if  Bi(t)  has  m eigenvalues  with  real  parts  < — 2i/  < 0 for  each 

te[o,T]. 


\Qi(t,s)\  < k3e-’*t—Vt  0 <s<t<T  (4.42) 

|02(*,s)|  < 0<^<s<r  (4.43) 

where  k3  is  a positive  constant  independent  of  e. 

Let  4>(C0)  be  the  fundamental  matrix  such  that 


Ht,s)  = A(t)^t,s)  4>(S,S)  = /2„X2„ 


45 


Then,  |4>(f,.s)|  is  bounded  by 


\${t,s)\  < k4eki l‘-l 


(4.44) 


for  some  constant  k4  and  k5  [29]. 

Using  the  transition  matrix  <!>(£,  s)  of  the  slow  system,  we  rewrite  (4.3S)  and 
(4.40)  into  a convenient  form.  Since 


rTN+tNy 


—\r  — r1 

ai {Tn  + £ 7*)  — »i(T/v)+  al(s)ds 

JTn 


^n(^,0)a1(0)  + $12(f^,0)a2(0) 
rfN 


+ Jq  ($n (Tn,s)1{  + $i2(Tn,s)1^  ds 

rTN+cN1- 

+ 6t!  (s)ds 

J1N 

, — at  — 1 i* 

P(Tn  + e 7*)  = P(Tn)  + - [ ( B/3  + m*)ds, 

t JfN 

equation  (4.38)  can  be  rewritten  by 

«i(0) 

$n(TV,  0)a!(0)  + $i2(2at,  0)q2(0)  + x^-(7)v)7* 


L fUv+dV/^O/  s -O  ,rf\  u , 

tN  JTn  \xn{s)  — xN[TH))ds 


+ eP1SM0{  0) 


+eP2SM  \P{TN)  + i J^N+t  7 ( B/3  + m^ds'j  - P2  J*"  <$>(TN,s)lmds 
[■Tn+(.N1* 

-P-i  {Aa  + l*)ds  + 0(e~ATN+cN-y')/c\ 

Jtn  j 


(4.45) 
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Ilcnce,  combining  (4.45)  with  (4.40)  gives 


where 


ai(0) 

h\(a,(3, 7*) 

L(T„,e) 

«2(0) 

= 

h2(a,(3, 7*) 

7* 

h3(<*,P,Y) 

(4.46) 


L(TN,e)  = 


L31 

L32 


10  0 
0)  4>12(X)v,0)  x^(Tn) 

L31  L32  C3 

— Ci$u(Ttf,  0)  + C2$2\(Tn,  0) 

= Cf  <&i2(Tn,  0)  + Cj ^22(?)v,  0) 


h\(a,P,  7*) 

0 

_ Ma»0.7*)  _ 

— 7^  If^+C  ^ (%n(S)  ~ ^n(^v))^s  _ 

(4.47) 


r*N 

+ePlSM/3(0)  — P2  $(TN,s)lmds 
Jo 

( _ 1 rTN+tN  7*  N 

+eP2SM  ( f3(TN)  + -J_  ( D(3  + m*)ds 


rTtt+cN -y* 

-P2  I (Aa  + l*)ds  + 0(e-'(T"+t"’y')/l)  (4.48) 
Jtn 


/i3(a,0,7*)  = -H*  + e[c?  C2r]  {PXSM  + P2SM)/3(fN) 

T 

- [cl  cj]  jf  "(Pi  + P2)$(TN,s)rds.  (4.49) 


Lemma  4.5.  There  exists  e3  > 0 such  that  for  all  0 < e < e3,  i)  the  (2n  + 1)  x 
(2n  + l)  matrix  L(T /v,  e)  defined  in  (4.46)  has  a bounded  inverse  and  ii)  given  K > 0, 
the  2m  x 2m  matrix  4/(7  , e)  defined  in  (4.41)  has  a bounded  inverse  if  |7*|  < Kt. 
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Proof  Assumption  A. 3 implies  that  the  linear  homogeneous  variational  TPBVP 


■o 


o — Anx  o + A^Pq 


jo 


Po  — A2iXq  + A22Po 


o 


0 


(4.50) 

(4.51) 


with  the  trivial  boundary  conditions 


xo(0)-  0 (4  52) 

xo  (To)  + ig  {T0)Ti  = 0 (4.53) 

0(^0)  + HZp%(T0)  + (MTo)  + Hj£°(To)  + HTpp°{TQj)  r, 


= 0 


(4.54) 


has  only  the  trivial  solution,  where 


A{t)  = uy-uwU-1Uy. 


Since 


xo(To)  = $u(To,0)xg(To)  + $ia(7b,0)p?(ro) 
Po(To)  = *2i(To,0)x°(To)  + <1>22(To,0)p°(To), 
the  boundary  conditions  (4.52)  — (4.54)  are  equivalent  to 


xo  (0) 

0 

T(T0, 0) 

o) 

= 

0 

Ti 

0 

Suppose  that  L(T0,  0)  is  not  invertible.  Then,  any  non-trivial  solution  start- 
ing with  a non-zero  initial  condition  satisfying  (4.55)  is  a solution  of  the  TPBVP 
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(4.50)  - (4.54).  Hence,  we  proved  that  L(T0,  0)  is  invertible  by  contradiction.  Since 
all  the  entries  in  L(TN,e)  are  continuous  functions  of  e,  there  exists  an  e3  > 0 such 
that  L{TN,t)  has  a bounded  inverse.  This  completes  the  first  part  of  the  proof. 
From  (4.41),  A/( 7*,e)  can  be  written  by 


m(7v) 


A/ii(0)  0 

0 M22(TN  + e^T*) 


+ 0(e). 


Since  Mn( 0)  and  M22(T0)  are  invertible  from  Assumption  A. 5,  M( 7*,  e)  is  invertible 
for  sufficiently  small  e.  a 


4.3  Main  Result 

Let  C(I,7l  ) be  a Banach  space  of  continuous  functions  on  a compact  interval 
I with  values  in  TV'  endowed  with  a norm 


|*(0loo  = sup  |ar(£)| 
tel 

where  / = [0,  TN  + eN+1I{]  for  some  positive  I<  and  t.  We  define  a Banach  space 
W with  a norm  |io|e  by 

w = c(i,n2n)xc{i,n2m)xn 

IHe  — kl(0lco  + |w>2(0loo  + |^3| 


where 


xu  = 


wi(t) 

(0 


W3 


w^t)  ec(i,n2n)  w2{t)  £C{i,n2m)  w3eiz. 


For  an  n x m matrix  A(t),  ||;4(f)||  denotes  suptg/  |/i(<)|. 
Consider  a closed  subset  B of  W defined  by 
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B = B,  x B2  x B3 

where 

e,  = {x(()eC(/,K2”)  | |x(<)U<A'e} 


8,=  W‘)€C(a")  I |x(()U<Ae} 

B3  = {x  £ 71  | |x|  < Ke). 


Now,  we  construct  a nonlinear  operator  $ acting  on  B such  that  to  solve  the 
TPBVP  (4.36)  and  (4.37)  with  the  boundary  condition  (4.38)  - (4.40)  is  equivalent 
to  find  a fixed  point  a,  /3  and  7*  of  the  nonlinear  operator  'L.  Thus, 


$ 


*i  (0 
^2(0 
r3 


: B ->  C(7, 7Z2n)  x C(7, 7Vm)  x 71. 


For  7r  (E  Bi,  £ 6 Bi  and  p E £>3,  we  define  T such  that 


r : B ->  7l2n  x 7l2m  x 71 


Tx 

ri(7T,  <^,  /?) 

hi{n,£,p) 

r = 

r2 

= 

r2(7T,£,/>) 

= L(TN,e)~1 

r3 

T3(t r,£,p) 

M7r>£>/>) 

Since  7y(T)v,e)  has  a bounded  inverse  (Lemma  4.5),  we  can  obtain  initial  conditions 
Fi  and  T2  of  the  slow  variable  and  a final  time  T3  from  a map  (4.56). 

1 herefore,  a trajectory  of  the  slow  variable  can  be  determined  by  ^i(^) 

*i(f)  : B^C(I,R2n) 


®i(0  = ’Mf*  *»£./>) 
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= *M) 


r2 


+ [ <b{t,s)r(TT,£,p,e)ds.  (4.57) 

JO 


Since  M(p,e ) is  invertible,  a boundary  condition  for  the  fast  variable  can  be 
determined  from  (4.39). 

A : B -*  n2m 


A = 

Ax 

Ai(tt,£,p) 

a2 

_ A2(tt,£,p)  _ 

r3,e)-x  (QM'i(O)  + (/  - + eNT3)) 


— A/(T3,  e) 


-l 


^ rfr"''*  j 

Since  a boundary  condition  for  the  fast  variable  is  determined,  we  obtain 
trajectory  for  the  fast  variable  by  defining  a map  \P2(0  as  follows. 

MO  • 

^2  (<)  = ^2{t,TT^,p) 

0i(<,O)A1 
®2  + eArr3)A2 


+ 0(e-^+tNr3)/£).  (4.58) 


1 

+- 

e 


Jo  Oi(t,  s)ml(Tr,£,  p,  e)ds 


(4.59) 


/iw+c^rs  02(4>  £>  P,  e)ds 

Theorem  4.0-  Suppose  that  Assumptions  A.l  - A. 5 hold.  There  exist  e*  and  5* 

such  that  when  e € (0,e*]  and  K(0)|,  \w*(0)\  < 6*,  the  solution  of  the  TPBVP 

(4.17)  and  (4.18)  with  the  boundary  condition  (4.21)  - (4.25)  exists  and  the  solution 

y*,  w*  and  7*  are  0(e). 
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Proof.  It  suffices  to  prove  that  the  nonlinear  operator  T has  a fixed  point  in  a 
subset  of  B.  Our  proof  consists  of  two  steps.  First,  we  will  show  that  |T(t,  7r,£,  p)|e  < 
Kt.  Next,  we  will  show 


TV  — 7T 

i-i 

p-p 


for  all 


■ 

■ 

TV 

TV 

i 

? 

i 

1 

1 

p 

G B. 


By  a Banach  space  fixed  point  theorem  [26],  this  implies  that  the  operator  T has  a 
unique  fixed  point  in  a subset  of  B. 

Assume  that  T defined  in  Lemma  4.3  is  sufficiently  large  such  that  T > TV  + 
eN+1K . Throughout  the  proof,  for  each  *,  * = 0,1,2,-*-,  c,  denotes  a constant  which 
is  independent  of  K and  u(e)  denotes  a continuous  function  such  that  lim^o  f (c)  = 0. 

Recall  that  xj)  = e~vt^  + e-v(TN-iN+'K-t) /e  jn  Lemma  4.2.  Therefore,  for  any 

t € I 

f xpdt  < kee2uKlN  for  some  k. 

Jo 

As  a preliminary  step,  from  (4.44) 

\\$(t,s)\\  < k4ek^TN+tN+lKl 

For  any  finite  K > 0,  the  choice  of  e4  = l/K  gives  ||<P(f,s)||  < Co  for  all  0 < e < e4 
and  t,s  E I.  Hence,  using  Lemma  4.4  and  |7r|,  |f  |,  |p|  < Ke  gives 


1/  $(t,s)l*(Tv,£,p,e)di 
\Jo 
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rTN+iN+'K 


||$(<,s)|||r(7r,£,p,e)|c/s 

rTN+i^+'K 


<-l 

r1 

— L kl  (c  + W + (e  + ‘W’XM  + |£|)  + + 60ip / e)\p\)  ds 

J o 

< (cj  + c2e2uKtN(80K  + 6)  + v(e))  e.  (4.60) 

Using  (4.42)  and  Lemma  4.4, 

7 1/ 

1 I U 


- 7 \Jo  k^e~‘/(t~a)/t  (e  + (e  + 5o)(k|  + \t\  + |p|))  ds 


k\k$ 


< —(1+3  K80  + v(t))t 


v 


< (ci  + c2e2^N(80K  + S)  + v(c))  e. 


(4.61) 


Similarly, 
1 I f* 


7 |XN+^+iA'02^,5)m2(7r,^^C)^'S  - (ci + c2e2i/A£N(^oA"  + <5)  + t;(e))  e.  (4.62) 

In  equations  (4.48)  and  (4.49),  since  |l£(s)  - i%(TN)\  < k\s  - TN\  for  some  jfc, 


*N 


rTN+cNp  _ 

(x%(s)  - xN(TN))di 

J In 


< §e"|,|2 

k 


< V+2A'2  = u(e)e 


^Ts+tN  p 

Jtn 

IT" 

By  using  (4.60)  and  (4.63)  - (4.65),  one  can  estimate 


rTN+eNp 

L (£|£|  + m*)ds 

J 1 N 

(yl|7r|  + /*)ds 


< u(e)e 

< v(e)e. 


(4.63) 

(4.64) 

(4.65) 


hi 


< (ca  + c4e2l/KeN (8 0K  + 8)+  v(e))e. 


53 


Since  L(TN,t)  has  a bounded  inverse  (Lemma  4.5),  |r(7r,f,/j)|  is  bounded  by 


< \\L(Tn,  e)-1!! 


h\ 

M7 r,€,p) 


h3  (V,{,p) 


< (c5  + cee2l'Kc  ( S0K  + 6)  + u(e)j  e. 


(4.66) 


From  a map  (4.57),  using  (4.60)  and  (4.66)  gives 


7T,^,/E»)|oo  < ||$(<,0)|| 

ri(7T,^,/9) 

+ 1/  $(t,s)T(Tr,Z,p,e)ds 

_ r2(7T,^,p)  _ 

\J  0 

< (c7  + c8e2"Ae  (80I<  + 8)  + u(e))  e.  (4.67) 


‘From  (4.58),  |A(tt,  £,p)|  can  be  written  as 


|A(7T,£,/>)|  < ||M(r3,e)-1|||g^1(0,7r,ie,p)| 

+ ||M(r3,e)-1|||(/-g)^1(fiV  + e;vr3,7r,e,p)| 
+^||A/(r3.£)-‘|| 

Since  |^(t)|  is  a known  function  which  is  bounded  in  the  compact  interval 
[0,  Tn  + eN+1I{],  using  Lemma  4.2  and  (4.67)  we  obtain 


rTN+tNr3(iT^,p,i)  Q 

qN(s)ds 

J±N 


+ 0(e 


-vTn/c\ 


|A(7t,^,p)|  < (09  + cme2i'KcN (80K  + 8)  + u(e))  e.  (4.68) 


Finally,  from  (4.61),  (4.62)  and  (4.68)  the  map  |'fr2(^,  tt,  £,  p)\oo  in  (4.59)  is 
bounded  by 


hMoo  < 


0l(f,O)A1(7T,^,p) 

®2 (t,  Tn  + eArr3)A2(7r,  £,  p ) 


OO 
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7/o©i(M)™;(7r,£,/>,e)<k 
[ 7 frN+c»T3  ®2{t,s)mm2(TT,Z,p,e)ds 

J oo 

< (cn  + C\2e2l/ht  (60 K + 8)  + u(e)j  e.  (4.C9) 

Let  Ci3  = c5  + c7  + C\ i and  c14  = eg  -|-  c8  + c12.  Then 

W,*,  t,  P)  le  < (ci3  + clAevKtN  (80I<  + 8)  + v{c))  e 
We  choose  sufficiently  small  e5  and  8 such  that 

C13  + ci4e2i/Ae  (80K  + 6)  + u(e)  < K for  all  0 < e < e5  0 < 80  < 8 


Consequently, 

l'L(C7T,£,/>)|e  < Ke 

This  completes  the  first  part  of  the  proof. 

Throughout  the  remainder  of  the  proof,  if  h is  a function  of  w,  £ and  /?,  then 
h and  h represent  h(ir,^,p)  and  /i(7f,  <f,  p),  respectively.  Also  d{  for  i = 0, 1,2,  ••• 
denotes  a constant  uniformly  bounded  in  e. 

As  a basic  step,  by  using  Lemma  4.4 

I [' -l*)ds 
'Jo 

,,  rTN+cN+1  K 

- ll$(L-s)||  Jo  k2( \tt  - 7r|  + |£  - £|  + |p  - p\)ds 

< k2c'0( \tc  - Trloo  + I i - eu  + e2^! p - p\) 


< k2c[( Iff  - *!«,  + I i - + | p-  p\)  (4.70) 

- I / ©i(*,s)(mr  - rnl*)ds 

C'J°  oo 

- l\l  - *1  + I i - ll  + 1/5-  p|)ds 
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< k2c'2 (|tt  - + |f  - ll^  + I P-  p|). 


(4.71) 


Similarly, 


7 \JfN+(Nt3  02(<’  ^ ^4(1*  - t|oo  + \i  - lu  + |p-  p~|).  (4.72) 

By  using  the  basic  estimates  (4.70)  - (4.71),  we  can  obtain 

|r(*,£p)-r(7r,£,p)|  < fc2C;(|7r  - Trloo  + - ||oo  + |p  - p|)  (4.73) 

l^i ~ *,Lp)\oo  < he's (|tt  - + | l-IU  + |p-p|)  (4.74) 

|A(tt,|,p)- A(ir,(f,p)|  < fc2c;(|7T-7r|00  + ||-e|oo  + |p-p|).(4.75) 

The  details  of  this  computation  are  long  and  tedious  but  straight-forward  and  so 
they  will  be  omitted  here. 

1^2  — 44 |oo  is  obtained  as  follows.  Since  from  (4.43)  |02(t,s)|  is  0(1)  for  all 
t,s  € I,  |02(M)|  is  0(1 /e)  from  (4.42).  Hence, 


||02(t,  TN  + e^fs)  - 02(f,  TN  + eNf3)|| 

< -j-  |f  - f I < A:2  — (1 7T  - 7f|oo  + | £ - (floo  + |p  - p|). 


(4.76) 


Since  m2  is  0(e), 
1 rfN+tN  r3 


r * n~t£  i 3 

eJTN+t„ra  H02(^4IIK*|ds  < 4| f 3 - f3 


< hc'10(\*  - TfU  + I i - lu  + |p  - p|)-(4.77) 
By  using  (4.71),  (4.72)  and  (4.75)  - (4.77),  |T2  - is  bounded  by 
| T2  4,2|oo  < ll0i(TO)|||Ai  — Ai|  H — |/  Oi(t,  s)(m\*  — rrii*)ds 


T||02(T T]\f  + eNV3)  - Q2(t,TN  + ejVf 


3 Pv2 
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+ ||02(<,Tyv  + eAf3)|||A2  - A2|  + - / . ||02(f,s)|||m2*|(/.s 

e Jtn+iN r3 

+^l/wf,0!(<’s)(n5j'“rtj')<isL 

< fc2c'n(|7r  - xloo  + |£  - fU  + Ip  - p|).  (4.78) 

Hence,  I'P  — ^|e  is  bounded  by 

I®  “ ®|e  < ^2ci2(|7T  - 7r|oo  + |£  “ f |oo  + \f>  ~ p\) 
where  c'12  = 04  + c'5  + c'u. 

Now,  we  choose  k2  such  that  k2dx2  < 1.  Finally,  let 

e*  = min(e!,  e2,  €3,64,65^! (fc2))  6*  = min(5,  Xi(k2)) 

Then  for  e < e*  and  <50  5-  £*,  the  nonlinear  operator  ^ has  a unique  fixed  point  in  a 
subset  of  B.  This  completes  the  proof.  ■ 

4.4  Concluding  Remarks 

In  the  proof  of  the  main  result,  Assumptions  A.l  - A. 5 are  made  on  the  lower 
order  problem.  The  fact  that  these  lower  order  assumptions  allow  the  determination 
of  the  asymptotic  solution  is  the  utility  and  power  of  the  singular  perturbation 
approach.  We  require  that  the  size  of  the  zeroth  order  boundary  layer  jumps  80  to 
be  sufficiently  small.  But,  this  is  inevitable  in  a general  nonlinear  context.  Similar 
results  can  be  found  in  Freedman’s  work  [22]  which  considered  a nonlinear  free  end 
point  optimal  control  problem  with  a fixed  final  time. 

In  proving  our  main  result,  we  need  three  key  assumptions;  1)  sufficient  smooth- 
ness of  our  system,  2)  the  reduced  solution  satisfies  Jacobi’s  sufficient  condition  for  a 
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local  minimum  and  3)  an  eigenvalue  condition  for  the  fast  subsystem.  But  these  as- 
sumptions are  not  strong  in  our  opinion.  The  smoothness  assumption  is  required  to 
deiive  an  asymptotic  series  solution  up  to  any  degree  of  accuracy.  Jacobi’s  sufficient 
condition  guarantees  a well-posedness  of  the  reduced  problem.  Thus,  without  Ja- 
cobi’s condition,  the  auxiliary  optimization  problem  is  meaningless.  The  eigenvalue 
condition  for  the  fast  subproblem  is  inevitable  for  the  application  of  the  singular 
perturbation  method  to  a control  problem  since  we  can  not  neglect  the  higher  order 
fast  dynamics  if  the  fast  dynamics  are  unstable.  But,  we  only  assume  that  the  sys- 
tem matrix  of  the  fast  system  (in  our  notation  Uw ) is  not  a stable  matrix  as  a time 
varying  system  but  a pointwise  stable  matrix. 


CHAPTER  V 

NONLINEAR  SYSTEMS  WITH  MULTIPLE  PARAMETERS 


In  this  chapter,  we  consider  a nonlinear  TPBVP  which  has  multiple  parame- 
ters. In  this  case,  we  can  construct  an  asymptotic  solution  by  a similar  method  as 
described  in  Chapter  III  in  several  time  scales. 

First  we  consider  a nonlinear  TPBVP  with  two  small  parameters 


X = 

fix,y,z) 

(5.1) 

ci! / = 

y(x,y,z) 

(5.2) 

cie2i  = 

h(x,y,z) 

(5.3) 

with  boundary  conditions 

xi(0)  = xlo 

Xi(T)  = XiT 

2/i(0)  = 2/i0 

to 

II 

o 

(5.4) 

o 

II 

O 

*7 

z2(T)  = 0 

and  a transversality  condition 

H(x,  y, 

N 

if 

S 

II 

o 

(5.5) 

where  Xj,  y i and  z\  are  n,  ni\  and  m2  dimensional  state  variables,  x2,  y2  and  z2 
are  n,  m\  and  m2  dimensional  costate  variables  corresponding  to  x1?  yi  and  Z\, 
H{xiy->z)  is  a scalar  function  and  e,  for  i = 1,2  is  a small  positive  parameter. 
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The  above  TPBVP  (5.1)  - (5.5)  comes  from  application  of  Pontryagin’s  min- 
imum principle  to  a final  time  unspecified  optimal  control  problem  for  a nonlinear 
system  with  two  small  parameters.  O’Malley  developed  an  asymptotic  expansion 
method  in  an  initial  value  problem  for  a nonlinear  system  with  two  small  parameters 
[30],  Also,  Khalil  investigated  stability  of  the  full  system  from  stability  of  the  re- 
duced system  in  a large  scaled  nonlinear  system  with  multiple  small  parameters  [31] 
[32].  In  this  section,  we  will  extend  the  asymptotic  expansion  method  developed  in 
Chapter  III  to  a nonlinear  singularly  perturbed  TPBVP  with  multiple  parameters. 

In  (5.1)  - (5.5),  by  setting  ej  and  t2  to  zero,  we  obtain  the  reduced  problem 


XT  = 

f(xri  Vri  zr) 

(5.6) 

0 = 

Uti  zt) 

(5.7) 

0 = 

Vri  zr) 

(5.8) 

*lr(0)  = Xl0 

3-lr  (Tr)  = XiT 

H(^xr , 

Vri  zr)\i-j'r  = 0. 

The  objective  in  this  section  is  to  seek  an  asymptotic  series  solution  based  on 
the  reduced  solution.  In  order  to  further  proceed,  we  need  some  assumptions  similar 
to  the  earlier  assumptions  mentioned  in  Chapter  III. 

Assumptions. 

HJ-i  f i 9i  h and  II  are  infinitely  differentiable  with  respect  to  all  their  argu- 
ments. 


ITT  In  equation  (5.8),  zT  — cf)(xr,yr)  is  a continuously  differentiable  solution 
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satisfying 


0 = h(xr,yr,<f>(xr,yr)) 


Let  yr  — i/>(xT)  be  a continuously  differentiable  solution  such  that 


0 = g(xr,il>(xr),<l>(xr,il>(xr))) 


The  reduced  problem  (5.6)  - (5.8)  has  a unique  and  infinitely  differentiable  solution 
on  the  interval  0 < t < Tr. 

B.3.  The  reduced  solution  satisfies  the  sufficient  conditions  for  a local  minimum. 
B.4.  There  exists  a 2m!  x 2 continuously  differentiable  matrix  M(t)  for  all 
0 < t < Tr  such  that 


M (tydyr  (■£»•)  Utj  ^{Xf,  yr^A'I (i)  — 


B{t)  0 

0 - DT{t ) 


Moreover,  Afn(0)  and  M22(Tt)  are  nonsingular  and  all  the  real  parts  of  the  eigen- 
values of  B{t ) are  negative  for  all  0 < t < Tr. 

B.5.  There  exists  a 2m2  x 2m2  continuously  differentiable  matrix  N(t)  for  all 
0 < t < Tr  such  that 


N l{t)hZr(xr,yr,zr)N(t) 


C(t)  0 
0 -CT{t) 


Moreover,  iVu(0)  and  N22(Tt)  are  nonsingular  and  all  the  real  parts  of  the  eigenval- 
ues of  C(t)  are  negative  for  all  0 < t < Tr. 

To  this  end,  we  will  obtain  an  asymptotic  double  expansion  in  and  e2  in  live 
different  time  scales.  We  will  define  four  stretched  time  scales  by 


t t 

Ti  — — T2  = <7 1 

el  ^1^2 


T-t  T-t 

(72  = 

Cl  £1^2 
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We  want  to  seek  an  asymptotic  series  solution  of  the  form 

x(t,e  1,62)  = x°(t,  ci,  62)  + £i  (xLl  (rj,  ej,  e2)  + xRl(cri,  ci,  £2)) 

+eie2  (xL2(t2,  Ci,  e2)  + xR2(a2,  ex,  e2))  (5-9) 

y(<>CiiC2)  = + (yLl(Ti,e1,e2)  + yRl((rueU€2)) 

+e2  (j/^(t2,  ex,  e2)  -f-  yR2(cr2,  cx,  e2))  (5.10) 

z(t,ti,e2)  = z°(t,ei,e2)  + (^zLl  (Ti,ei,t2)  + zRl(cri,ei,e2)^ 

+ [zL2(t2,  Ci,  e2)  + zR'2(a2 , ei,e2))  (5-11) 

Here,  all  the  terms  on  the  right  hand  side  of  (5.9)  - (5.11)  have  an  asymptotic  double 
power  series  expansions  in  and  e2  as  both  e.\  and  e2  go  to  zero.  For  example, 

x°(t,eue2)  = x°k(t)44 

j,k>0 

t = Ew  4. 

j.fc>0 

In  this  case,  j/qo)  zoo  an<^  ^00  will  satisfy  the  reduced  problem  (5.6)  - (5.8). 
In  the  outer  region,  by  substituting  x°(i,  ex,  e2),  y°(<,e1,e2)  and  z°(t,el,e2)  into 
(5.1)  - (5.5)  and  equating  coefficients  of  equal  powers  of  for  j + k > 0,  we  have 


T° 

Xjk 

~ fx(x00,y00,z00)xjk  + fy(x00,y00,  z00)yjk 

+ fz{x 00,  2/00 > Z00 )Zjk  d"  ^1  ,jk(t) 

(5.12) 

0 

— 9x(x00,  y00,  z00)xjk  + gy(x00,y00,z00)yjk 

~\~9z{x00i  2/001  ZOo)Zjk  "b  ^2 ,jk(t) 

(5.13) 

0 hx(x00,  y00,  z00)x°k  -f-  hy(xQ0,  y°0,  ZQQ)y^k 


+ hz(x 00,  y00,  ^oo)2jl  + ^3,jfc(0 

0 = ^(xoo)2/oo»2roo)xifc(^oo)  + + R°jk 
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(5-14) 
(5.15) 

where  i?°fc  for  1 < i < 5 are  known  successively  in  terms  of  x°n,  y°n,  and  Tmn 
with  m + n < j + k. 

With  the  aid  of  Assumptions  B.4  and  B.5,  z°k  can  be  obtained  in  terms  of  x°k 
and  yfk  from  (5.14)  and  yfk  can  be  obtained  in  terms  of  x°k  from  (5.13).  In  order 
for  an  asymptotic  solution  to  satisfy  the  boundary  condition  (5.4),  we  require  that 
for  j + k > 0, 


0 

0 


<(°)  + <M(°)  + 4t1.»-,(o) 

(°)  + KiA  + Kik 


(5.16) 

(5.17) 


where  and  R% jk  are  successively  known  values  and  quantities  with  a negative 
index  are  defined  to  be  zero. 

In  the  first  left  boundary  layer^  scale),  we  will  substitute  (5.10)  - (5.11)  into 
(5.1)  - (5.3)  assuming  that  the  terms  in  r2,  <J\  and  cr2  scale  are  negligible.  Then  the 
first  boundary  layer  terms  must  satisfy 
For  j = k = 0, 

dxLl 

dTl  ~ f(xoo(0),yoo(0)  + 2/oo  (ri))  ^00(0)  + •2:oo(ri)) 


dyoo 

d,T\ 


(5.18) 


/(•Too(0)>  2/oo(0)i  zoo(0)) 

5(^o(0),y?o(0)  + yoo(n),4(0)  + 4>(ti)) 
-5(^S)(0),yg,(0)^°(0)) 


(5.19) 
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0 — Mxoo(0))  2/oo(0)  + Voo  (ri)i  2oo(0)  + zoo  (ri)) 

-M^(0),yo°o(0),4(0)) 

For  j + k > 0, 
dxLl 

(jT^  = /v(xoo(0)>  J/oo(0)  + Voo  (Ti)>  zoo(0)  + zoo  (ri))2/jfc  (Ti) 

+ /z(xoo(0)>  J/oo(0)  + Voo  (Ti)>  ~oo(0)  + zoo  (Ti ))zfk  (ri) 

+Ri!jk(r1) 

dyLl 

jTi  = ^(^oolO))  2/oo(0)  + yoo(ri)?  2oo(0)  + •2:oo(Ti))2/ji1(Ti) 

+ /z(xoo(0)>  J/oo(0)  + Voo  (Ti  )j  zoo ( 0 ) + zoo  (ri  ))zjk  (Ti) 

+R2Jk(Tl) 

0 = ^^(^00(0)?  J/oo(0)  + Voo  (ri))  zoo(0)  + zoo  (ri))j/jfc  (ri) 
+/i2(x°(0),y°(0)  + J/oo  (Ti)>  zoo(0)  + zoo(Ti))z:k(Ti) 
+ R3Jk(Tl) 

with  the  initial  condition 

’jt(O)  = yu  - J,°o(0) 

= -y?J0)  - y^_t(0) 


(5.20) 


(5.21) 


(5.22) 


(5.23) 


where  Ri}jk(ri)  f°r  1 < * < 3 are  successively  known  functions  and  decay  to  zero  as 
Ti  goes  to  infinity. 
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From  (5.20)  or  (5.23),  zj'/j  can  be  obtained  in  terms  of  y^  algebraically.  Hence 
(5.19)  and  (5.22)  is  a differential  equation  containing  only  yfy . The  initial  condition 
of  the  fast  costate  ( 0 ) must  be  so  chosen  that  yfy  decays  exponentially  to  zero 
as  T\  goes  to  infinity.  Then,  since  the  right  sides  of  (5.18)  and  (5.21)  are  known, 
can  be  obtained  by  integration.  The  exponential  stability  in  the  first  left  boundary 
layer  is  guaranteed  if  Assumption  B.4  holds  and  |t/f^(0)|  is  sufficiently  small. 

Similarly,  in  the  second  left  boundary  layer(r2  scale),  the  second  left  boundary 
layer  terms  must  satisfy 


For  j = k = 0, 


<fopo 
dr 2 

dy0p 

dr2 

dzpp 

dr2 


/(xoo(0)>  2/oo(0)  + 2/oo  (0)?  zoo(0)  + zoo  (0)  + zoo  (r2)) 
-/(*&(0),y£(0)  + yti(0),z°(o)  + 2oo  (0)) 
^(^oolO))  Z/oo(0)  + Voo  (0)>  2oo(0)  + 2oo  (0)  + zoo  (r2)) 
-^o(0),yo°o(0)  + 2/oo(0)^o°o(0)  + 2oo(0)) 
M*8>(0),  Voo(°)  + Voo  (°)>  zoo(° ) + zoo  (°)  + zoo  (r2)) 
— /i(xoo(0),  j/Oo(0)  + y0Q  (0),  2oo(0)  + 2oo(^))- 


(5.24) 


(5.25) 


(5.26) 


For  j + k > 0, 


dxfk 

dr2 


f z(xoo(0),  J/oq(0)  + J/oo  (0))  2oo  + zoo  (0)  + zaS(T2))zjk(T*) 
+<^(r2)  (5.27) 

^(■rOo(0)i  !/oo(0)  + Voo  (0),  Zqq  + Zqo  (0)  + zoo(T2))zjk(T2) 


+#2,yr2) 


(5.28) 
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dzfk 

dr2 


^z(;roo(0))  2/oo(0)  + J/oo(0)i2oo  ~b  2oo(0)  + 200(T2))~jfc(T2) 
+^:k(r2)  (5.29) 


with  the  initial  condition 


*£(°)  = *..-*?„(  o)-*&(o) 

4;»(0)  = -2?;1(0)-zf;,(0) 

where  R^k(r2)  for  1 < i < 3 are  successively  known  functions  and  decay  to  zero 
as  r2  goes  to  infinity.  By  Assumption  B.5,  we  can  choose  z%*k( 0)  such  that  zfk(r2) 
decays  exponentially  to  zero  as  t2  goes  to  infinity  if  |z22fc(0)|  is  sufficiently  small. 

By  the  same  manner,  the  first  and  the  second  right  boundary  layer  expansion 
can  be  obtained.  Moreover,  extension  of  this  asymptotic  expansion  method  to  a 
nonlinear  TPBVP 


ele2 • • • tkVk 


y(®>yi»---»yn)  i<^<« 


with  n small  positive  parameters  ek  for  1 < k < n will  be  possible  by  a straight- 


forward manner. 


CHAPTER  VI 

APPLICATION  OF  SINGULAR  PERTURBATION 
TO  BTT  MISSILE  GUIDANCE 


6.1  Problem  Statement 


Preliminaries. 

The  proportional  navigation  (PN)  guidance  law  has  been  widely  used  for  a 
missile  guidance  during  the  fast  several  decades  because  of  its  effectiveness  and 
simplicity  of  implementation.  In  the  PN  guidance  law,  the  acceleration  command 
is  generated  by 

ac  — A nav^c^ 

where  ac  is  the  acceleration  command,  vc  is  a closing  velocity,  a is  a line  of  sight 
(LOS)  angle.  The  navigation  constant  Knav  is  determined  experimentally.  It  was 
shown  that  the  PN  guidance  law  was  optimal  under  certain  conditions  [25].  The  PN 
guidance  law,  however,  is  not  generally  effective  against  maneuvering  targets  [7]. 

In  order  to  obtain  better  performance,  many  authors  have  attempted  to  derive 
a guidance  law  by  using  modern  control  theory.  Actually,  the  guidance  problem  for 
an  air-to-air  missile  may  be  formulated  as  an  optimal  control  problem.  Since  the 
missile  motion  equation  is  nonlinear,  the  resulting  nonlinear  TPBVP  is  difficult  to 
solve  by  on-line  computation. 

In  order  to  circumvent  this  computational  burden,  Fiske  developed  an  optimal 
linear  (OL)  guidance  law  by  using  a simplified  missile  motion  equation  modelled  as 
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double  integrator  in  the  inertial  axis  [33].  If  a missile  is  modelled  as  a point  mass  in 
the  inertial  axis,  the  missile  motion  equation  can  be  casted  into  the  following  form. 


where  x,  vx  and  ax  are  the  relative  distance,  relative  velocity  between  the  missile 
and  a target  and  an  acceleration  command  in  the  inertial  X-axis. 

The  OL  guidance  law  was  obtained  by  minimizing  the  following  quadratic  per- 
formance index. 


where  Tgo  is  computed  by  an  appropriate  time-to-go  estimation  algorithm.  The  in- 
eitial  X-axis  acceleration  command  can  be  easily  obtained  in  a closed  form.  By  the 
same  method,  the  inertial  Y-axis  and  Z-axis  acceleration  command  can  be  com- 
puted. From  these  inertial  acceleration  commands,  an  acceleration  command  in 
the  body  axis  can  be  obtained  by  the  inertial  to  the  body  axis  coordinate  transfor- 
mation. The  desired  roll  angle  command  is  computed  such  that  the  body  Y-axis 
acceleration  command  is  zero.  Since  the  body  X-axis  acceleration  is  nearly  constant 
in  our  concerned  BTT  missile,  the  computed  body  X-axis  acceleration  command 
should  be  neglected  in  Fiske’s  OL  guidance  law.  The  performance  of  the  optimal 
linear  guidance  law  is  satisfactory.  However,  when  the  initial  range  is  short  (within 
4000  ft),  the  performance  of  the  OL  guidance  law  is  not  satisfactory.  This  may  be 
due  to  the  neglected  missile  body  X-axis  acceleration. 

Recently,  an  attempt  has  been  made  to  derive  a nearly  optimal  guidance  law 
from  a nonlinear  missile  motion  equation  by  the  singular  perturbation  method. 


x 


v 


v. 
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Sridar  derived  a nearly  optimal  guidance  law  which  is  approximated  up  to  the  zeroth 
order  [7].  But,  Sridar  observed  that  the  solution  significantly  deviates  from  the 
optimal  as  time-to-go  decreases. 

In  order  to  obtain  a better  approximation,  Visser  suggested  a first  order  cor- 
rection method  by  using  the  matched  asymptotic  expansion  method  [9].  But,  it 
was  observed  that  the  nearly  optimal  guidance  law  derived  based  on  the  singular 
perturbation  method  also  failed  to  give  a good  approximation  during  the  terminal 
phase  in  which  the  time  scale  decomposition  is  not  valid. 

Problem  Statement,. 

We  define  the  reference  axis  as  the  wind  axis  at  each  sampling  time.  Then,  the 
missile  motion  equations  in  the  reference  axis  will  be 


dt 

dyw 

dt 

dzw 

dt 

d~/ 

dt 

di/> 

dt 

d(j> 

dt 

with  the  initial  condition 


V cos  7 cos  0 — Vt 


txvu 


V cos  7 sin  -0  - Vtyw 
Vsini  - Vtzw 

I 

cos  0— 


sin  0 


V cos  7 


Pu 


(6.1) 

(6.2) 

(6.3) 

(6.4) 

(6.5) 

(6.6) 


x™(0)  — xqu,  Vw(0)  = y0w  zw{  0)  = zqw 


7(0)  = 0 0(0)  = 0 0(0)  = 0 


where 
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xw,yw,zw  : the  relative  distance 

Vtxw,  Vtyun  Vtzw  : the  target  velocity  (assumed  to  be  constant) 

V : the  missile  average  velocity 

The  flight  path, heading  and  roll  angle 
Lw  :the  vertical  lift 
pw  :the  roll  rate 

\Lyj\  ^ Lwmax,  \pw\  < Pwmax- 

In  this  chapter,  our  aim  will  be 

1.  to  derive  a nearly  optimal  guidance  law  with  the  first  order  correction  terms 
using  the  method  described  in  Chapter  V. 

2.  to  increase  the  region  in  which  the  time  scaling  decomposition  is  valid  by 
adjusting  control  weighting. 

6.2  Singularly  Perturbed  TPBVP 

This  section  contains  a procedure  to  transform  (6.1)  - (6.6)  into  a standard 
form  of  a singularly  perturbed  system.  In  order  to  formulate  the  problem  in  the 
form  for  singular  perturbations,  we  normalize  the  variables  and  time  scale  as  follows. 
Let 

R = (*L  + vL  + zD* 

X = Xw/R  y = yw/R  z = Zw/R 

^ tr  - V txw/R  Vty  = Vtyw/R  Vtz  — Vtzw/R 
L ~ Lw/ Lwmax  p = pw/pwmax  t'  = Vt/R. 
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Using  the  above  definitions,  equations  (6.1)  - (6.6)  become 


x = cos  7 cos  ip  — Vtx 

(6.7) 

y = cos  7 sin  ip  — Vty 

(6.8) 

z = sin  7 - Vtz 

(6.9) 

/ii7  = cos  p L 

(6.10) 

Piip  = sin  <p——— 
cos  7 

(6.11) 

Pip2(p  - P 

(6.12) 

with  the  initial  condition 


where 


a:(0)  Xqw/R  — xq  y(0)  — yow/di  — Vo  ^(0)  — ZqwJ  R - zo 
7(0)  = o V>(0)  = 0 </>(0)  = 0 


x = 


Pi  = 


dx 

dt' 


V2 


RLr 


P2  = 


PmaxV 


The  controls  L and  p will  be  chosen  to  minimize  a quadratic  performance  index 
with  a cheap  control  cost. 


J = [ + \k\V  + \k\kl?)it 

subject  to  the  final  state  and  control  constraints  for  an  unspecified  T . 

*)  <T)  = y(T)  = z(T)  = 0 
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»)  1^1  <1  \P\  < 1 


where  0 < k2  <C  1. 

The  Hamiltonian  for  this  problem  is 


H = Ax(cos7  cos - Vtx)  + Aj/(cos7  sin^  — Vty)  + A2(sin7  — Vtz) 

+A7  cos  (f>  L + Xj,  sin  <f>— - b A,*  p 

cos  7 

+ 1 + 2 + 2 ^ ^2  P2 

where  for  example,  Ax  denotes  the  costate  corresponding  to  x. 

Now,  we  define  new  controls  and  new  costates  as  follows. 


L — L p — fc\  Ic2  p \x  — A x 

Ay  = Ay  A2  = Az  A7  = — — 

k\ 

\ — ht  \ — ^ 

~ ki  X<t>  ~ klk*' 

Then,  from  the  Euler-Lagrange  equation,  the  necessary  conditions  for  optimal- 
ity will  be 


H 


x 

y 


Ax(cos7  cos  xj>  — Vtx)  + Ay  (cos  7 sin0  — Vty)  + A2(sin7 

+ A7  cos  <f>  L + A^  sin  <f>— b A ^ p + \l 2 + \p2  + 1 

cos  7 2 2 


Vu) 


0 

cos  7 cos  xk  — Vtx 
cos  7 sin  0 — Vty 
sin  7 - Vtz 


(6.13) 

(6.14) 

(6.15) 


z 


(6.16) 
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Cl  7 

= 

cos  p L 

(6.17) 

ci  P 

- 

sin  pL 
cos  7 

(6.18) 

Cl  C2  4> 

= 

P 

(6.19) 

^ x 

o 

II 

II 

3) 

— -<C 

(6.20) 

Ci  A7 

= 

dH 

d~f 

(6.21) 

= 

dH 

dip 

(6.22) 

Ci  c2  A ^ 

dH 

dp 

(6.23) 

L 

- 

arg  min  H 
IU<*i 

(6.24) 

P 

= 

arg  min  H 

\p\<k\k2 

(6.25) 

where 

-if 

II 

e2  = n2  h 

x(0)  = x0 

O 

II 

o 

si 

2(0)  = Zq 

7(0)  = 0 

o 

II 

o 

O 

II 

O 

Remark 

6.1. 

Actually,  the  accuracy  of  an 

asymptotic  solution  depends 

on  the 

size  of  the  singular  perturbation  parameter  e. 

But  in  an  actual  system  to  be  con- 

trolled,  it  is  difficult  to  identify  the  exact  value  of  the  singular  perturbation  param- 
eter. Usually,  a designer  can  estimate  a singular  perturbation  parameter  roughly 
based  on  his  insight  into  a dynamic  system  or  rough  estimates  of  the  time  derivatives 
of  the  states  as  in  deriving  (6.7)  - (6.12). 
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Remark  6.2.  In  equations  (6.7)  - (6.12),  the  singular  perturbation  parameter 
and  /r2  rnay  not  be  valid  during  the  whole  interval.  For  example,  /q  becomes 
bigger  as  the  time-to-go  goes  to  zero  (the  relative  distance  is  shorter).  In  order 
to  compensate  rough  estimates  of  ji\  and  /j, 2,  we  introduced  a performance  index 
with  a cheap  control  cost.  In  the  TPBVP  (6.13)  — (6.25),  the  singular  perturbation 
parameter  t\  and  e2  are  functions  of  /ii,  //2,  &i,  By  choosing  small  k\  and  fc2,  the 
singular  perturbation  parameter  and  e2  can  be  arbitrarily  small. 

Remark  6-3.  Introducing  a small  control  weighting  in  a performance  index  can 
be  found  in  earlier  papers.  O’Malley  and  Jameson  utilized  this  scheme  for  a singular 
linear  quadratic  optimal  control  problem  [34]  [35].  In  order  to  seek  an  optimal 
solution  in  a singular  optimal  control  problem,  they  regularized  the  problem  by 
introducing  £ in  a control  cost  and  obtained  an  asymptotic  solution  by  the  singular 
perturbation  method.  Ardema  also  investigated  the  behaviour  of  the  solution  of 
a nonlinear  singular  control  problem  on  the  singular  arc  by  regularizing  the  initial 
problem  [36]. 


6.3  Derivation  of  the  Guidance  Raw 

In  this  section,  we  will  derive  an  asymptotic  solution  of  the  nonlinear  TPBVP 
(6.13)  - (6.25)  in  a closed  form.  In  obtaining  the  asymptotic  solution,  the  most 
difficult  part  of  the  solution  procedure  is  to  solve  the  zeroth  order  boundary  layer 
equations  whose  solution  must  be  asymptotic  to  the  zeroth  order  outer  solution. 
Once  the  zeroth  order  solution  is  obtained,  the  first  order  correction  terms  can  be 
made  based  on  a procedure  described  in  Chapter  III. 

A linearization  approach  [37] [38]  for  solving  the  nonlinear  zeroth  order  boundary 
layer  equations  was  suggested.  But,  in  this  case,  when  initial  conditions  of  fast  states 
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are  far  away  from  the  zeroth  order  outer  solution,  there  is  no  guarantee  such  that 
the  initial  conditions  are  in  the  region  of  attraction.  When  i)  nonlinear  boundary 
layer  equation  is  scalar,  ii)  the  final  time  is  unspecified  and  iii)  the  Hamiltonian 
is  not  an  explicit  function  of  time,  a closed  form  solution  can  be  easily  obtained 
from  the  fact  that  the  Hamiltonian  is  identically  zero.  In  this  section,  we  will  use  a 
coordinate  transformation  such  that  the  2 dimensional  boundary  layer  equations  in 
the  first  left  boundary  layer  can  be  treated  as  a scalar  boundary  layer  equation. 

Outer  Solution. 

By  setting  ea  and  e2  to  zero  in  the  outer  expansion,  (6.14)  - (6.23)  become 


•O 

xoo 

cos 7°  cos V>oo  - Vtx  *oo(0)  = x0  x°(T00)  = 0 

(6.26) 

Voo 

= 

cos  7o°0  sin  V>oo  - Vty  Voo(0)  = y0  y°o(Too)  = 0 

(6.27) 

Zoo 

= 

sin7oo  ~ Vtz  ^oo (0)  = zo  z°0(T00)  = 0 

(6.28) 

0 

= 

cos  (PoqLqq 

(6.29) 

0 

= 

smPoo  00  o 

cos  7qo 

(6.30) 

0 

= 

Poo 

(6.31) 

\o 

/vx00 

= 

\0  \0  rv 

Ay00  ~ Az00  ~ U 

(6.32) 

0 

= 

^roo  s'n  Too  cos  'Poo  d"  Ay00  sin  700  sin  i/’oo  ~ ^^oo  cos  Too 

\0  ■ ±o  sin  Too  To 

'Voo  Sln  Poo  cog2  ^Oq  Loo 

(6.33) 

0 

- 

^x00  c°s  Too  sin  ''Poo  - tfoo  cos  7oo  cos  'Poo 

(6.34) 

0 

= 

^xoo  Sln  Poo  Loo  ~ \oo  cos  <P°oo  °° o ■ 

cos  7oq 

(6.35) 
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From  (6.34),  t/’oo  's  constant.  From  (6.29)  and  (6.30),  L°0  = 0-  Then  from 
(6.33),  7^0  is  also  constant. 

Hence,  from  (6.26)  - (6.28),  we  can  compute  T00,  J00  and  V>oo>  be., 


Too 

-C  + VC2  + AB 
A 

(6.36) 

7o°o  (t) 

_ . _1  VtzToo  — Zq 

— tan  - 

sJ(VtxToo  — %o)2  + (f^yToo  — Do)2 

(6.37) 

Oo 

- tan-1  VtyT°°  ~ y° 
VtxToo  — Xq 

(6.38) 

where 


A = l-(vi  + V,l  + Vtl) 

B = *0  + Vo  + xl 

C — Vtxx0  + vtyy0  -f  Vtzz0. 


We  can  obtain  further  simplification  by  redefining  the  roll  origin,  that  is  we 
perform  a rotation  about  the  X-axis  by  an  angle  A <f>.  Then, 


" 

" 

X 

X 

y 

= T 

y 

z 

z 

T = 


where 

1 0 0 

0 cos  A cj>  sin  A 0 

0 — sin  A (j)  cos  A 0 
In  the  new  coordinate  system,  7,  0 and  0 are  defined  as  before  with  respect  to 
the  new  axis.  Since  the  functional  form  of  the  equations  is  unchanged,  we  have 


0oo  = tan 


- f_-l  VtyToo  - Vo 


KxToo  — Xo 
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By  a choice  of  A <f>,  we  can  force  xpQQ 


= 0 as  follows. 


VtyToo  — iio 


(cos  A<pVty  + sin  A(j>Vtz)T0o  - (cos  A<py0  + sin  A <pz0) 


= 0 


if  A(f>  is  chosen  to  satisfy 


A<p  - — tan 


-1  Vo  — VtyTp 

zo  — VtzTo 


In  the  new  axis, 


7oo(0  = tan 


-i  VjzToo  — zo 
VtxToo  %o 


(6.39) 


7(0)  = 0(0)  = 0 ^(0)  = -A  <f>. 


In  the  following,  we  choose  the  new  variables  by  x,  y,  z , 7,  ip  and  <p  in  place  of 
x,  y,  z,  7,  ip  and  for  simplicity. 

Also,  the  H = 0 equation  becomes 


1 + A °oo(cos  Too  cos  0oo  - Vte)  + Ay00(cos  7°0  sin  V»S>  “ Vf») 
+^S)o(sin7oo  - Vfz)  = 0 (6.40) 


From  (6.33),  (6.34)  and  (6.40), 


Ao  _ ~ cos  Jq0 

1 -(Vtx  cos  Jqq  + Vtz  sin  7&) 

A°oo  s 0 


(6.41) 

(6.42) 


\0  _ \0  1 o 

^200  Ax00  tan  700. 


o 


(6.43) 


First  Left  Boundary  Layer  Solution. 

Since  jg0  in  (6.39)  does  not  satisfy  the  initial  condition  7(0)  in  general,  boundary 
layer  analysis  is  required  to  eliminate  this  discrepancy. 
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Since  rp°0,  <f>w,  and  p°0 


. w are  identically 

following  set  of  the  zeroth  order  necessary  conditions  in  the  first  boundary 

r L, 

^oo  — cos  Too 


rlrLl 

= cos(7to  + 7oLo  ) cos  ^0- 

^Voo  / o , L,\-  I 

-j^T  = cos(7oo  + 7oo  ) sin  0OO 


e?Zf 


A 

'00 


dri 


(1tx 

^7oc i 
d,Ti 

iM 


M 


sin(7oo  + Too  ) _ sin  7^> 


d^xoo  _ d\yQ0  — dAfo0  ^ 


drx 

COS  </>QQ  Z/QQ 


= sin  L™ 

cos(7o*+7,k) 


0 


Poo 


d\Ll 

dTl  ~ "^a-oo  sin(7oo  T Too  ) cos  V’oo  ~ Azoo  cos(7oo  T Too* ) 

Koo  sin  <f>00  sin(7o°o  + Tod  ) 2,  o°° 

cos2(7o°0  + 7oo‘ ) 


d\Ll 

= Ax00  cos(7°  + 7£j  ) sin  0, 


’oo 


0 


A700  s’n  ^oo  -^oo  A^do  cos  ^00  , q°°  ZifT 

cos(7o°0  + 7od ) 


Too1  (0)  = 7(0)-7o°o(0) 


-7o°o(0) 


V’oo  (0)  = V’(O)  - 0oo(°)  = 0. 

From  (6.47),  we  have 


(6.44) 

(6.45) 

(6.46) 

(6.47) 

(6.48) 

(6.49) 

(6.50) 


(6.51) 

(6.52) 

(6.53) 


A*do  — 0 Aydo  = o A(-n  = 0 


z00 


(6.54) 


Since  0oo  (0)  — 0,  no  boundary  layer  transition  can  be  expected.  This  fact 
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motivates  the  coordinate  transformation  described  earlier  in  this  section. 
^od(ri)  = 0 in  (6.49)  gives  V’oo(Ti)  = 0.  Hence  by  using  (6.54)  and 

V’od(ri)  ~ ^od  (Ti)  = 0? 


the  Hamiltonian  (6.13)  in  the  first  boundary  layer  will  be 


Hoo  = A*oo(cos(7od  + 7oo)  ~ Vtx)  + Af00(sin(70L01  + 7oo)  ~ Vtz) 
+^700-^00  + 1 + - (^oo)  = 0. 

Since  H°0  = H& , combining  (6.40)  and  (6.55)  gives 

A700^od  + 2 {Loo)  = Ax0o(cos7oo  ~ cos  (7^  + 7°)) 

+A?oo(sin7°o  - sin(7oo‘  + 7°o )) 

-2 


• 2 7oo 
sin 


From  (6.24), 


where 


1 ~ (Vtx  cos  7q0  + \/(3  sin  7^)  2 


/ 


1 if  /(x)  > 1 

sat  (/(x))=  /(x)  if  -1</(x)<1 

-1  if  f(x)  < -1. 

Hence  from  (6.56)  and  (6.57),  the  unsaturated  control  Lqq  satisfies 


V 


T L\  T/  • 700 

T'oo  = A sin  — 


I<  = 


2 


Setting 


(6.55) 


(6.56) 


(6.57) 


(6.58) 


where 
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If  the  control  is  saturated  (i.e.,  \K  sin  > k\), 


Lqo  = ~ki  s’Sn(7oo  )• 


(6.59) 


Since  Lqq  is  known  by  (6.58)  or  (6.59),  A^q0  can  be  computed  from  (6.56) 

al,  \ sin2  + \{LqqY 


7OO 


L 


(6.60) 


oo 


Since  4)  = 0,  from  (6.53)  we  obtain 


Second  Left  Boundary  Layer  Solution. 

In  the  second  left  boundary  layer,  application  of  (5.24)  - (5.26)  using 


xP 


o 

00 


l — \°  — \0  _ \L\  _ lO  _ /Li  _ \0  _ \L\  _ TO  _ rv 
TOO  — 'SOO  ~ AipOO  — 'Voo  — 'Poo  ~ YOO  ~ A<\> 00  — A^00  ~ ''00  — u 


gives 


dxoo 

dr2 

d^xoo 

dr2 

d^od 

dr2 

dip  op 
dr2 

dX  700 
dr 2 

d^T^o  0 

dr2 

d<Poo 

dr2 


dy0o 

dr2 


= 0 


d~od 

dr7 


= 0 


d\Pl 


yOQ 


dr7 


= 0 


dxpi 


zOO 


dTo 


= 0 


cos  4o2  (4o  (0)  + 4o)  - Loo  (0) 


sin  4o2 


Lqq (0)  + 
!(7oo(0)  + 


0 

0 


Poo 


dXjlp 

dTo 


= A^oo(O)  sin4o(-^oo(0)  + 4o) 


(6.61) 

(6.62) 

(6.63) 

(6.64) 

(6.65) 

(6.66) 

(6.67) 


(6.68) 
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with  the  initial  condition 

^oo  (0)  = -A  (f). 

From  (6.65)  and  (6.66), 

A->Oo(T2)  = \oo(0)  A^0(t2)  = ^V'Oo(O)  = 0. 

Hence,  from  the  Hamiltonian  in  the  second  left  boundary  layer 
Loo  = arg  min  II ^ 

= -*isat  cos  - Lfo  (6.69) 

Poo  = arg  min  II ^ 

Poo 

= -W2Sat(i^)-  (6-7°) 

Since  the  Hamiltonian  is  constant,  from  II ^ = H^j  we  obtain 

2 (Poo  ) + A^oqPoo  = ~A^q0  cos  <^oq  (Fqq  + Fqo)  + A^q0Lqo  . 

-j((ife+ife)2-(iS)2)  (6.71) 

Using  (6.69)  and  (6.70),  from  (6.71) 

Poo  = -fci^2sat(y— )sign  {$*)  (6.72) 

M /^2 

where 


= ^oocos^o2  (/&>  + Z^)  - A^qZqq  + i ((Z&  + Zqq)2  - (Zqo)2)  . 


We  remark  that  if  the  zeroth  order  control  is  implemented  by  a feedback  scheme, 
all  the  zeroth  order  state  variables  appeared  in  (6.69)  and  (6.72)  can  be  replaced  by 
the  exact  values  since  all  the  state  variables  are  measurable. 
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First  Order  Correction. 

Up  to  this  point,  the  zeroth  order  solutions  in  the  outer  region  and  the  first  left 
boundary  layer  region  are  derived.  These  solutions  are  derived  under  the  assumption 
that  the  slow  variables  are  fixed  while  the  fast  variables  reach  to  their  steady  state 
values.  But  this  assumption  may  not  hold  during  the  terminal  phase.  Since  the 
singular  perturbation  parameter  e2  is  a function  of  the  maximum  allowable  limits 
of  control  variables  and  the  missile  average  velocity,  e2  is  nearly  independent  of  a 
missile-target  geometry.  Hence,  we  will  derive  the  first  order  correction  terms  in 
in  order  to  increase  the  region  in  which  the  singularly  perturbed  solution  is  valid. 
As  can  be  seen  in  (6.56),  (6.57)  and  (6.69),  the  control  variables  are  functions  of 
only  the  slow  costate  variables  if  all  the  state  variables  are  measurable.  Hence,  it 
suffices  to  derive  the  first  order  correction  terms  in  the  slow  costate  variables. 

To  this  end,  for  simplicity,  we  will  define  xa  = Xq0-\-  e\X*0  where  the  superscript 
a is  0 or  Ll.  In  the  outer  region,  (6.14)  - (6.23)  can  be  expressed  up  to  the  first  order 
in  ex  in  the  following  form  if  is  sufficiently  small. 


x = cos  70  cos  xj>°  — Vtx 

y = cos  70  sin  4>°  - Vty 

•o  . ~ 

2 = sin  7 — Vtz 


0 

0 = sin  <jj 

0 = p° 


cos  4>° L° 
L° 


o 


cos  70 


(6.73) 

(6.74) 

(6.75) 

(6.76) 

(6.77) 

(6.78) 


(6.79) 
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0 = sin  70  cos  + A°  sin  70  sin  0°  — A°  cos  70 
* - L° 

-Ajsin^sin-i0— (6.80) 

0 = A®  cos  7°  - A®  cos  7°  cos  (6.81) 

0 = A°  sin  4>°L°  - Ag  cos  j°  ■ ^ (6.82) 

cos  7U  v 7 

The  TPBVP  (6.73)  - (6.82)  has  the  same  functional  form  as  the  TPBVP  in 

(6.26)  -(6.35).  The  initial  conditions,  for  example,  i°(0)  = 2^(0)  + €ixfo(0)  can 

be  determined  from  (6.44)  by  integration.  Since  the  right  hand  side  of  (6.44)  is  a 

known  function,  we  obtain  by  integration 


4(0)  - Jo  [cos (7^  + 7°)  cos (^oo1  + ipoo)  ~ cos 7oo cos V’oo] 

roo  r 

= / [cos(7oo  + 7oo)  - cos  7®  drx 

•/ 0 


_ [°  cos  (7^  + Too)  - cos 


= r 

J-T& 


'Li 

Joo 


00 


o 


Sm  7 sat  7sat  COS  7oo 

Mgnfrft) 


+ f cos  7^  - cos  75a<  + /0°)  (6.83) 


where  7 sat  is  defined  by 


if  l^oo1(0)|<l  7«t  = 0 


if  Lqq  (0)|  > k\  K sin 


7 sat  700 


if  Tqo(0)|<—  ki  A' sin 


^ . Isat  7oO 


—k\. 


Similarly, 


y?o(0) 

4(0) 


= 0 


_ - COS  7 aat  -f  1 - ~faat  sin  4 


^isign(70°0) 


+ J7  ( sin  7°  - sin 


• Isat  + 700 


. (6.84) 
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Since  x°(0),  y°{ 0)  and  3(0)  are  known  and  the  TPBVP  (6.73)  - (6.82)  have 
the  same  functional  form  as  those  of  the  TPBVP  (6.26)  - (6.35),  the  slow  costate 
variables  can  be  computed  by  the  same  method  described  in  the  derivation  of  the 
zeroth  order  outer  solution. 

We  remark  that  the  right  boundary  layer  terms  are  identically  zero  in  this 
particular  problem. 


6.4  Implementation 

Derivation  of  the  Body  Axis  Roll  Rate  Command. 

Since  the  equations  of  motion  (6.1)  - (6.6)  are  defined  in  the  wind  axis,  all 
the  initial  conditions  and  the  target  velocity  must  be  evaluated  in  the  wind  axis. 
Furthermore,  since  pw  in  (6.6)  is  defined  as  the  wind  axis  roll  rate,  it  must  be 
converted  to  the  body  axis  roll  rate  in  order  to  be  used  as  the  roll  rate  command  to 
the  autopilot.  The  wind  axis  roll  rate  can  be  expressed  in  terms  of  the  body  axis 
roll  rate  pb  as  follows 


Pw  = cos  /?  cos  Q pb  + sin  (3  qb  + cos  /3  sin  arb-  a sin  (3  (6.85) 


where  pb,  qb  and  rb  denote  the  body  axis  roll  rate,  pitch  rate  and  yaw.  a and  (3  are 
the  angle  of  attack  and  the  side  slip,  respectively. 

Since  (3  and  rb  are  negligibly  small  in  a bank  to  turn  missile,  pb  can  be  expressed 
as  follows 


Pb 


Pw 

cos  a 


(6.86) 


Determination  of  the  Control  Weighting  Aq  and  k2. 

The  validity  of  the  singular  perturbation  method  depends  on  the  insignificance 
of  the  higher  order  terms.  Since  ea  = as  the  final  time  is  approached,  the 
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missile  approaches  the  target,  R approaches  zero  and  tx  increases.  When  e,  increases, 
the  neglected  higher  order  terms  become  significant,  causing  the  approximation  to  no 
longer  be  valid.  This  is  the  reason  for  the  failure  of  the  usual  singular  perturbation 
method  in  the  terminal  phase.  To  avoid  this  problem,  kx  can  be  adjusted  to  cancel 
the  effect  of  R going  to  zero.  A simple  solution  is  to  pick  kx  so  that  ex  is  constant. 
A second  cause  of  failure  for  the  singular  perturbation  method  is  the  size  of  the 
coefficients  multiplying  ex  and  62  in  the  expansion.  If  these  coefficients  are  very 
small,  the  solution  will  remain  valid  even  as  ex  and  e2  get  large.  For  example,  if 
7(0)  ~ 7oo>  then  the  higher  order  coefficients  are  all  small  so  that  larger  values  for 
and  e2  do  not  destroy  the  approximation.  Due  to  modeling  errors  and  unmodelled 
target  motion,  7(0)  is  not  well  approximated  by  7$,  and  the  fixed  weight  singular 
perturbation  method  fails. 


CHAPTER  VII 
SIMULATION  RESULTS 


7.1  Deterministic  Simulation  Results 

The  performance  of  the  guidance  law  developed  in  the  previous  chapter  will  be 
demonstrated  by  using  6 DOF  bank-to-turn  missile  simulation  program  developed 
at  the  University  of  Florida.  Table  7.1  shows  the  missile  characteristics  used  in  the 
simulation  program. 


Table  7.1  Summary  of  Missile  Characteristics 


Initial  velocity  (Mach) 

0.9 

Boost  thrust  (lb) 

4712 

Burn  out  time  (sec) 

2.6 

Initial  weight  (lb) 

165 

Final  weight  (lb) 

115 

Maximum  vertical  lift  (g) 

100 

Maximum  roll  rate  (deg/sec) 

500 

The  simulated  target  is  a three  dimensional  maneuvering  target  with  the  max- 
imum acceleration  9 g.  Once  the  range  from  the  missile  to  the  target  is  less  than 
6000  ft,  the  target  initiates  an  instantaneous  9 g evasive  maneuver  in  a plane  de- 
termined by  the  target  roll  angle  which  is  determined  by  an  input  parameter.  The 
Rockwell  autopilot  which  uses  a dithered  signal  in  the  yaw  channel  was  used  in  the 
simulation  program.  There  is  also  a 0.4  second  delay  between  missile  launch  and 
autopilot  action. 


85 


86 


I he  miss  distances  of  the  singularly  perturbed  guidance  law  with  the  control 
weight  adjusting  scheme  (WASP),  with  a constant  control  weight  (WUSP)  and  the 
optimal  linear  guidance  law  (OL)  are  shown  in  Figure  7.1. 

From  Figure  7.1,  it  can  be  seen  that  the  weight  adjusted  singularly  perturbed 
guidance  law  developed  in  Chapter  VI  performs  well  under  the  various  engagement 
scenarios  in  spite  of  the  modelling  errors  such  as  the  target  maneuvering  and  ne- 
glected autopilot  missile  air  frame  dynamics. 


Figure  7.1:  Miss  distances 


In  order  to  illustrate  the  performance  of  the  missile  with  the  WASP  control 
foi  an  initial  range  of  5000  feet  and  initial  gimbal  elevation  and  azimuth  angles  of 
0 and  40  , is  shown  the  commanded  vertical  lift,  roll  rate  and  the  corresponding 
trajectory  are  shown  in  Figure  7.2  - Figure  7.5. 


Roll  Rate(deg/sec)  Lift(g) 


Figure  7.2:  Vertical  Lift  Command 


Figure  7.3:  Roll  Rate  Command 


Z position(feet)  Y positioned) 


Figure  7.4:  X-Y  Plane  Trajectory 


Figure  7.5:  X-Z  Plane  Trajectory 
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From  Figure  7.2,  the  vertical  lift  command  is  saturated  during  the  end  game. 
This  comes  from  the  fact  that  the  LOS  angle  to  the  intercept  point  must  be  zero  at 
the  final  time  to  satisfy  the  final  conditions  (the  relative  distances  at  the  final  time 
are  zero).  Hence,  when  the  time-to-go  is  extremely  small,  the  commanded  vertical 
lift  will  be  saturated  even  if  the  LOS  angle  is  small. 

In  order  to  further  compare  the  singularly  perturbed  guidance  law  with  perfor- 
mance index  adjustment  with  the  optimal  linear  guidance  law,  we  will  examine  some 
of  the  simulation  results  in  more  detail.  The  test  trajectory  chosen  for  illustration 
has  an  initial  range  of  10,000  feet  and  initial  gimbal  angles  of  0°  and  40°  for  elevation 
and  azimuth.  The  behavior  of  the  missile  wind  axis  line  of  sight  angle  (LOS)  to  the 
predicted  intercept  point  is  shown  in  Figure  7.6  for  the  weight  adjusted  performance 
index  (WASP)  and  the  constant  weight  performance  index  (WUSP). 

Both  schemes  keep  the  LOS  fairly  small,  indicating  good  tracking,  after  the 
initial  transient.  However,  the  constant  weighting  control  law  fails  to  track  in  the 
end  game  due  to  the  inability  to  track  the  unmodelled  target  maneuver.  In  all 
terminal  guidance  schemes,  the  time-to-go  is  a significant  parameter  which  must 
be  estimated.  The  estimate  given  by  the  range  divided  by  the  closing  velocity  has 
been  shown  by  Riggs  and  others  [39]  [40]  to  be  in  error  for  maneuvering  targets 
and  in  the  final  portion  of  the  trajectory.  Riggs  [39]  suggested  an  improvement, 
known  as  Riggs  algorithm,  which  takes  the  known  missile  acceleration  into  account. 
The  estimated  time-to-go  in  the  weight  adjusted  singular  perturbation  is  superior  to 
Rigg’s  approximation  as  shown  in  Figure  7.7  where  the  actual  time-to-go  is  compared 


with  the  two  estimates. 


Tg°(sec)  LOS(deg) 
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Figure  7.6:  LOS  Angle 


Figure  7.7:  Time-To-Go 
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Figure  7.9:  Side  Slip 
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In  Figure  7.7,  all  three  curves  were  generated  from  the  same  trajectory,  only 
the  method  for  determining  the  time-to-go  was  modified.  The  improvements  in  the 
time-to-go  estimates  may  be  another  reason  for  the  superior  performance  of  our  new 
guidance  law.  The  history  of  the  angle  of  attack  and  the  side  slip  for  an  initial  range 
of  5000  feet  are  shown  in  Figure  7.8  and  Figure  7.9.  From  Figure  7.9,  WASP  keeps 
the  side  slip  within  the  limit  (±5°)  for  a BTT  missile. 

7.2  Stochastic  Simulation  Results 

So  far,  we  assumed  that  the  target  velocity  is  known.  Actually,  the  only  avail- 
able measurements  with  respect  to  the  target  are  i)  the  range  measurement,  ii)  the 
range  rate  measurement,  iii)  the  elevation  angle  measurement  and  iv)  the  azimuth 
angle  measurement  when  an  active  seeker  is  used.  In  this  simulation,  the  extended 
Kalman  filter  was  used  to  estimate  the  target  velocity  from  the  above  four  avail- 
able measurements.  The  covariances  of  the  seeker  measurement  noise  used  in  the 
simulation  are  given  in  Table  7.2. 


Table  7.2  Covariance  of  Measurement  Noise 


Type  of 
Measurements 

Covariance 

Azimuth  and 
Elevation  Angle 
Range 
Range  Rate 

((1.5 /R)2  + 2.2252  x 10-10) /AT 

(9  + 10  ~16R4) /AT 
(9  + 4 x 10_22i?4)  /AT 

R:  the  relative  distance  (ft)  AT:  the  measurement  step  (sec) 


The  performance  was  determined  by  10  runs  for  each  initial  condition.  The 
mean  and  the  variance  of  the  final  miss  distances  with  the  initial  seeker  elevation 
angle  0.107°  and  azimuth  angle  40°  are  shown  in  Table  7.3 
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Table  7.3  Statistics  of  Miss-distances 


: Initial 

Guidance  law 

(ft) 

WASP 

WUSP 

OL 

10000 

Mean 

Var 

36.9 

749.2 

65.8 

604.8 

32.3 

572.9 

8000 

Mean 

Var 

24.0 

130.1 

44.8 

302.1 

24.5 

101.5 

7000 

Mean 

Var 

22.3 

126.5 

32.8 

107.1 

24.6 

159.2 

5000 

Mean 

Var 

30.8 

435.7 

30.9 

440.0 

34.8 

591.3 

4500 

Mean 

Var 

31.5 

309.4 

29.7 

186.4 

58.1 

285.8 

4000 

Mean 

Var 

41.9 

268.5 

44.3 

328.6 

75.3 

813.4 

3800 

Mean 

Var 

37.8 

295.8 

38.0 

339.4 

61.3 

701.8 

A typical  trajectory  and  command  history  of  an  initial  range  5000  feet  are 
shown  in  Figure  7.10  - Figure  7.13. 

The  average  performance  of  the  new  weight  adjusted  singular  perturbation 
control  law  is  superior  to  the  constant  weight  singular  perturbation  and  optimal 
linear  guidance  laws.  The  relatively  large  sample  variance  observed  for  the  miss- 
distances  indicates  more  runs  are  required  for  a more  thorough  evaluation. 


Z-axis  Y-axis 
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Figure  7.10:  X-Y  Plane  Trajectory 


Figure  7.11:  X-Z  Plane  Trajectory 
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Figure  7.12:  Vertical  Lift  Command 


Figure  7.13:  Roll  Rate  Command 


CHAPTER  VIII 
CONCLUSIONS 


In  this  dissertation,  the  singular  perturbation  method  was  considered  in  order 
to  analyze  a nonlinear  TPBVP.  Particular  emphasis  has  been  given  to  a nonlinear 
TPBVP  arising  in  a nonlinear  optimal  control  problem  with  free  final  time.  The 
dissertation  consists  of  two  main  parts. 

In  the  first  part,  an  explicit  method  to  construct  a uniformly  valid  asymptotic 
solution  of  the  full  problem  was  developed.  By  our  construction  method  based  on  a 
recursive  scheme,  the  solution  of  an  optimal  control  problem  can  be  approximated 
to  any  order.  With  the  assumptions  imposed  on  the  reduced  problem,  if  e and 
the  zeroth  order  boundary  layer  jump  S0  are  sufficiently  small,  the  solution  of  the 
full  problem  exists  and  approaches  the  asymptotic  solution  as  e goes  to  zero.  In 
the  proof  of  this  result,  one  of  the  crucial  assumptions  is  an  eigenvalue  criterion 
assumed  in  Assumption  A. 5.  However,  as  shown  in  Sacker’s  work  [41],  eigenvalues 
with  zero  real  parts  may  be  allowed  if  a form  of  conditional  stability  is  satisfied  so 
that  Assumption  A. 5 could  be  weakened. 

In  constructing  an  asymptotic  solution,  one  of  the  most  difficult  things  is  to 
obtain  the  zeroth  order  boundary  layer  solution.  However,  in  an  actual  application 
problem,  sometimes  it  may  be  sufficient  to  obtain  the  zeroth  order  asymptotic  solu- 
tion. In  this  case,  we  can  use  an  approximate  zeroth  order  boundary  layer  solution 
obtained  from  the  linearized  boundary  layer  problem  with  respect  to  the  zeroth 
order  outer  solution. 
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The  singular  perturbation  method  can  also  be  applied  to  a system  with  multiple 
time  scale  property  on  an  infinite  time  interval  [42].  One  of  the  interesting  questions 
in  this  area  is  that  if  a reduced  order  system  is  asymptotically  stable,  under  which 
conditions  a full  system  is  asymptotically  stable?  Recently,  this  question  was  partly 
answered  [43]  [44].  In  the  robustness  analysis  in  a control  problem  with  neglected 
higher  order  dynamics,  the  singular  perturbation  will  be  a powerful  tool. 

In  the  second  part,  the  singular  perturbation  method  was  applied  to  an  air- 
to-air  BTT  guidance  problem.  In  order  to  compensate  for  model  errors  and  target 
motion,  a control  weight  adjusting  scheme  is  used  to  derive  a guidance  law.  In 
both  the  deterministic  and  the  stochastic  cases,  this  scheme  gives  a considerable 
improvement  over  an  optimal  linear  guidance  law  when  the  initial  range  is  short. 
When  the  initial  range  is  long  (above  7000  ft),  this  scheme  also  works  well  while  the 
usual  singular  perturbed  guidance  law  with  constant  control  weight  does  not  work 
well.  The  computation  is  simple  enough  so  that  it  can  be  implemented  by  on  board 
computation. 

The  higher  gain  control  resulting  from  the  smaller  control  weighting  forces  the 
missile  to  quickly  reach  a coasting  arc  where  the  vertical  lift  command  is  zero.  If  we 
assume  that  i)  the  current  missile  position  is  on  the  coasting  arc,  ii)  the  missile  roll 
angle  can  be  changed  instantaneously  and  iii)  the  future  behavior  of  the  target  is 
unknown,  the  center  of  the  target  reachable  set  will  be  placed  on  the  center  of  the 
missile  reachable  set.  Thus,  the  idea  of  the  high  gain  control  implicitly  coincides 
with  the  idea  of  the  reachable  set  control  [45], 

The  performance  of  the  guidance  law  is  heavily  dependent  on  the  target  velocity. 
Since  the  covariance  of  the  measurement  noise  in  the  elevation  and  azimuth  angle  is 
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larger  as  the  missile  approaches  closer  to  the  target  due  to  glint  noise,  the  estimated 
value  of  the  target  velocity  is  not  good  during  the  terminal  phase.  This  leads  to  an 
increased  miss  distance.  In  order  to  further  reduce  the  miss  distance,  an  improved 
estimation  scheme  for  bearing  only  measurements  should  be  developed. 
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